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PREFACE 


When high-energy radiation is used to produce a change in matter— 
whether a physical change, as in the testing of materials for use in 
reactors, a chemical change, as in the polymerization of plastics, or a 
biological change, as in the treatment of cancer—the radiation must be 
measured and controlled. The measurement of radiation in situations 
such as these constitutes the subject of radiation dosimetry. 

Radiation dosimetry has always presented difficulties to those new to 
the field, and these difficulties have often stemmed as much from lack 
of understanding of the quantities being measured as from inadequacies 
in equipment or technique. The aim of this book is to present the special 
concepts of radiation dosimetry in a coherent manner and to show how 
they are related to such basic parameters of the radiation field as intensity 
and spectral distribution. The emphasis is on the physical principles 
of radiation measurements. When particular problems are mentioned 
for the sake of concreteness, they are mainly the traditional ones of 
biological and medical dosimetry, but in most cases the discussion applies 
equally well to the analogous problems in other fields. Although specific 
techniques and instruments are referred to occasionally by way of illus- 
tration, no attempt has been made to give a systematic or detailed 
description of radiation detectors. 

The book is intended primarily for the science graduate faced with 
the problem of determining radiation dose. Some familiarity with the 
elementary concepts of atomic physics and with simple differential and 
integral calculus is assumed; less familiar topics are discussed briefly 
as they are needed. Since an understanding of the various ways in which 
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radiations interact with matter is essential, this subject is treated in 
some detail in a separate chapter. 

A word about terminology is necessary. There is as yet no universal 
agreement on the names to be assigned to some of the physical quantities 
that occur in dosimetry. When possible I have followed the 1956 recom- 
mendations of the International Commission on Radiological Units and 
Measurements. Where these offer no guidance I have tried to adopt the 
usage of those most active in the field. Thus, for example, whereas the 
practice in other branches of physics would point to the use of the term 
“flux density” for the number of particles or photons passing through 
unit sphere per second, I have followed the usage common in neutron 
physics (and growing more common in dealing with other radiations) 
and used the shorter term ‘‘flux.’”? A number of such choices have been 
necessary, and some of them will almost certainly differ from what is 
eventually accepted. This is regrettable but unavoidable. 

Finally, I should like to express my thanks to those who read the 
manuscript and made helpful criticisms and suggestions: Drs. H. O. 
Wyckoff and R. S. Caswell of the National Bureau of Standards, 
Washington; Mr. F. H. Attix of the Naval Research Laboratory, Wash- 
ington; and Dr. W. R. Dixon and Mr. J. H. Aitken of the National 
Research Council, Ottawa. 

: G. N. WHYTE 


Ottawa, Canada 
April, 1959 
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CHAPTER 


1 


The radiation 


field 


When a beam of radiation traverses matter, three types of physical 
information may be of interest: 

(1) The spectral distribution of the radiation, 

(2) The intensity of the radiation (or, alternatively, the flux of parti- 
cles) at some point, 

(3) The amount of energy absorbed per unit mass at a point in the 
irradiated material. 

As will be seen further on, the ionizing effect of the radiation in air 
can sometimes be used to determine one or the other of the latter two 
quantities, and consequently a knowledge of the ionization in air under 
specified conditions may also be useful. 

Since much of the succeeding discussion will revolve around the prop- 
erties of radiation fields, it may be helpful to specify the meanings of the 
above concepts in some detail. 

Intensity and spectral distribution. For x-rays and 7-rays, the 
radiation field is completely specified if at every point we know the 
energy per second in the form of photons with energies between hy and 
hy + d(hv) crossing a unit area normal to the direction of propagation 
within the solid angle dQ in the direction defined by the angle variables 
6 and ¢ (see Fig. 1.1). This quantity may be written (hv, 6, ¢)d(hy) dQ, 
where I (hy, 6, ¢) may be termed the differential energy and angle dis- 
tribution of the radiation intensity. [Since J(hy, 0, ¢) will in general 
vary from point to point, it will also be a function of the position varia- 
bles z, y, 2.] 

A complete knowledge of this distribution is seldom needed. If it is 
summed over all directions (i.e., all values of @ and ¢) the result is 

I 


2 PRINCIPLES OF RADIATION DOSIMETRY 
I,(hv), the “spherical” energy (or spectral) distribution 
I,(hy) = f f I(hv, 6, 6) dQ (1.1) 
O“"¢ 


I,(hy) d(hv) can be thought of as the energy per second in the form of 
photons with energies between hy and hy + d(hy) passing through a 
sphere of unit cross-sectional area from any direction. This is what 
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Fig. 1.1. Illustrating the geometrical parameters involved in the definition of 
intensities. 
is detected by a device such as a scintillation counter which can measure 
photon energies without regard to direction. The rate of absorption of 
energy from the radiation field by matter is determined by I,(hyv). 

If the spectral distribution J,(hy) is further summed over all values of 
the photon energy hy, the result is J,, the spherical, or ‘‘omnidirectional,”’ 


intensity 
iF = { 1.0) d(hv) 


= J { I(hv, 8, 6) d9 d(hv) (1.2) 
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I, is the total energy per second of all photon energies passing through 
a sphere of unit cross-sectional area at the point of interest. [See Fig. 
1.2(a).] This is the quantity referred to as “intensity of radiation” in 
reference 1 and simply as “intensity” by most radiation workers. Be- 
cause of the recurring necessity of making clear distinctions between 
them, the terms “spherical intensity’ and “plane intensity’ (defined 
below) will be used throughout this book. In practice, the term “‘in- 
tensity” should be reserved for the former quantity. 





Fig. 1.2. (a) The spherical intensity is determined by the energy crossing unit 
sphere. (b) The plane intensity is determined by the energy crossing unit area. 


Where, as is often the case, the radiation is all travelling in more or 
less the same direction, it is sometimes useful to deal with ‘“‘plane’”’ or 
“unidirectional’’ intensities. The plane energy (or spectral) distribu- 
tion, I,(hv), is defined so that the energy per second in the form of pho- 
tons with energies between hy and hy + d(hv) crossing a unit area 
oriented in a particular direction is J,(hv) d(hy), 1.e., 


his J J, I(hv, 6, 6) cos ¥ do (1.3) 

6 
where y is the angle between the direction of propagation and the normal 
to the area. Unless the radiation is isotropic (in which case [,(hv) = 0), 


the value of I, (hv) will depend on the orientation chosen for the unit area. 
Summing over hy yields the plane intensity, J, 


I= f I,(hy) d(hv) 


= ff fre», 6, d) cos y dQ d(hy) (1.4) 
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I, is the total energy per second of all photon energies crossing a unit 
area at the point of interest. [See Fig. 1.2(b).] In the case of a parallel 
beam of radiation, the plane intensity, J,, referred to an area normal to 
the beam is equal to the spherical intensity, J,. Except for this particu- 
lar case, the spherical intensity always exceeds the plane intensity, as can 
be seen by comparing Eqs. 1.2 and 1.4. 

The time integral of the spherical intensity over a fixed interval is 
often of interest; it is termed the quantity of radiation. 

An example will illustrate the importance of the distinction between 
plane and spherical intensity. Consider a parallel beam of radiation, 
infinite in extent, striking a hypothetical medium in which the only in- 
teraction is scattering without loss of energy, i.e., the only effect of the 
medium is to change the direction of some of the incident rays. We shall 
take the direction of the incoming beam to correspond to 6 = 0, and con- 
fine our attention to the rays scattered through a particular angle, 6’, 
at any azimuthal angle, ¢. 

In Fig. 1.3 the rays that will be scattered through an angle & are shown 
coming into the medium from above and the scattered rays are shown 
emerging from below. If the incoming plane intensity is J,°, then the 
incoming spherical intensity, 7,°, is also equal to J,°. Let the scattered 
intensity in the direction 6, ¢ and in any interval d¢ of the azimuthal 
angle be J(@’) dd, and assume /(6’) independent of ¢. The plane intensity 
below the medium will be 


2r 2r 
I, = f 1(6") cos 6’ do = 1(6’) cos 6’ f de 
0 0 


= 2x1(6’) cos 6’ (1.5) 


But since there is no absorption, the same energy must cross any hori- 
zontal plane in the figure, and the plane intensity will be the same below 
the medium as above it. Therefore 


2x I(6’) cos 6’ = Ip° (1.6) 
The spherical intensity below the medium will be just 


2r 
Ee J 1(6") dé = 2nI(6") (1.7) 


Thus from Eqs. 1.6 and 1.7, 
I, = I,/cos 6’ = I,°/cos & 
= I,°/cos 6’ (1 8) 
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Whereas the plane intensity is the same above and below the medium, 
the spherical intensity is greater below. <A small ionization chamber 
would give a greater reading below the scattering medium than above. 

This situation is illustrated in Fig. 1.3. The two thick horizontal 
lines represent unit areas normal to the direction of the incoming beam; 
each is crossed by six rays, indicating that the plane intensities in the two 


Incoming beam 


AA 


Scattering medium 


Fig. 1.3. X-rays scattered through an angle @ in a non-absorbing medium. 





positions are equal. Similarly, the circles represent spheres of unit cross- 
sectional area. The one in the incoming beam is also crossed by six 
rays—the incoming plane and spherical intensities are equal—but the 
one in the outgoing beam is crossed by nine rays, indicating that here 
the spherical intensity is greater. 

Obviously the above arguments can be extended to the rays scattered 
at other angles, and hence to the entire scattered beam. The conclusion 
stands that in the absence of absorption the plane intensity will be the 
same at any depth, but the spherical intensity will increase with depth 
as an increasing fraction of the beam undergoes scattering. 

In any actual situation there will be some absorption along with the 
scattering, and the effects of this will be superimposed on the effects of 
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the scattering. Thus the plane intensity will decrease with increasing 
thickness of the medium, and the spherical intensity will either increase 
at first and then decrease, or simply decrease less steeply than the plane 
intensity. The effect of the scattering will still be to increase the spheri- 
cal intensity relative to the plane. 

Particle flux. A charged-particle field can also be specified in terms 
of intensities, using definitions similar to those set down above for x-rays. 
However, it is more common to consider the number of particles crossing 
unit area, rather than the energy. The flux, ¢, of charged particles at a 
point is defined as the number of particles per second crossing a sphere of 
unit cross-sectional area at that point.* (Flux is thus analogous to spheri- 
cal intensity.) When the particles have a distribution of energies, the 
energy spectrum is characterized by the differential flux, ¢(£), where 
¢(E) dE is the flux of particles with energies between E and E + dE. 


‘The total Hux in thia-casein piven by f $(E) dE, where the integral is 


taken over the range of energies present. 

The neutron field at a point is also commonly specified in terms of the 
flux. The neutron density, i.e., the number of neutrons per cm’, is fre- 
quently of interest as well. If the density of neutrons of velocity 
v cm/sec is n neutrons/cm®, the flux of such neutrons is simply 


@ = nv (1.9) 


Absorbed dose. The intensity and the related quantities defined 
above are characteristics of the radiation field. In most applications, 
however, what is of primary significance is not the energy in the beam 
but the energy absorbed per unit mass of irradiated material. A biologi- 
cal effect, for example, must depend on the amount of energy locally 
absorbed at the point in question, rather than on the amount of energy 
passing through this point. (It may, of course, depend on other factors 
as well, such as the rate at which energy is absorbed, or the photon energy 
of the radiation.) 

The quantity expressing the amount of energy locally absorbed is the 
absorbed dose, defined as follows (1): 


Absorbed dose of any ionizing radiation is the energy imparted to 
matter by ionizing particles per unit mass of irradiated material at the 
place of interest. 

The unit of absorbed dose is the rad. 1 rad is 100 erg/g. 


* This quantity might more properly be called flux density, but current usage favors 
flux. 
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The absorbed dose is not a property of the radiation field alone; it 
results from the interaction of the radiation with matter and therefore 
depends on the properties of the particular material involved. 
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The interactions 
of radiation 
with matter 


All radiation measurements depend on the interactions of the radiation 
with matter. The nature of these interactions must therefore form the 
basis of a discussion of the measurements themselves. This chapter con- 
tains an outline of the principal interactions of charged particles, x-rays, 
and neutrons with matter, with particular emphasis on the aspects that 
are of most concern in dosimetry. 


CHARGED PARTICLES 


Charged particles lose energy in matter chiefly through Coulomb inter- 
actions with atomic electrons. If the energy transferred to an electron 
is only enough to raise it to a higher energy level in the atom, the process 
is called excitation; if the electron is given enough energy to separate it 
completely from the atom, the process is called ionization. The two 
processes are closely associated, and together they constitute energy loss 
by collision. At energies below several times the rest energy of the mov- 
ing particle excitation and ionization account for the major part of the 
energy loss in all materials. 

Some of the electrons ejected in ionization processes have large enough 
energies to produce further ionization themselves; such electrons are 
called 6-rays. 

A number of other phenomena also occur when charged particles tra- 
verse matter. Interactions with the Coulomb fields of atoms, and par- 
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ticularly of atomic nuclei, result in changes in the direction of motion of 
the particle. If no kinetic energy is transferred to another form, the 
process is called elastic scattering; such scattering is of minor importance 
for heavy particles but of great importance for electrons. Scattering 
through a large angle entails a large acceleration of the charged particle. 
This in turn may result in the emission of a quantum of electromagnetic 
radiation, known as bremsstrahlung. 

If the speed of the particle through a medium exceeds the speed of 
light in the medium, radiation in the visible region known as Cerenkov 
radiation is emitted. 

Collision loss and stopping power. According to the classical 
point of view, a moving particle loses energy to an electron by imparting 
to it an impulse proportional to the strength of the Coulomb force and 
to the time during which this force acts. Since the momentum acquired 
by the electron is proportional to the time during which the interaction 
takes place, it is inversely proportional to the velocity, v, of the moving 
particle. The energy acquired by the electron, and hence the energy 
lost by the particle, must therefore be proportional to 1/v*. Thus on 
the classical picture the energy loss per unit length of path (specific 
energy loss) should be proportional to the electron density in the 
medium and inversely proportional to the square of the velocity of the 
particle. 

Because the electrons in matter are bound to atoms, the energy trans- 
fers do not occur out to infinite distances from the track of the moving 
particles. The farther an electron is from the track of the particle, the 
smaller the impulse it can receive, and hence the smaller the energy that 
can be transferred to it. If that energy is just less than the amount re- 
quired to raise a K-electron to a higher energy level, no energy will be 
lost to K-electrons at that distance. Somewhat farther away losses to 
L-electrons will become impossible, and so on. The more tightly bound 
the atomic electrons (i.e., the higher the atomic number), the shorter 
will be these “cut-off” distances, and the smaller the rate of energy loss. 
Consequently the rate of energy loss should show some dependence on 
atomic number, being less at high atomic numbers. 

At velocities approaching the velocity of light the 1/v? dependence is 
modified by a relativistic effect. The relativistic contraction of the 
electric field of the moving particle makes possible energy losses at 
greater distances, and in consequence the rate of energy loss increases 
slowly at very high energies. 

The crude arguments presented above serve only to bring out some of 
the general features of energy loss by collision. The quantum mechanical 
treatment has been given by Bethe (1, 2, 3). For heavy particles (i.e., 
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particles other than electrons) the specific energy loss, or stopping power, 
—dE/dz, is given by 


dE 4wz7e4 a Z [m Qmv? 6 4 ) (2.1) 
es — | In —————  — 8? — - | erg/cm 
dz met "AL I — 6?) 2) °° 


where ze is the charge of the moving particle in esu, e is the charge on 
an electron, m is the rest mass of the electron in grams, v is the velocity 
of the particle in cm/sec, No is Avogadro’s number, p is the density of 
the stopping material, Z is the atomic number of the stopping material, 
A is the atomic weight, and 8 is equal to v/c, where c is the velocity of 
light. The quantity J is the mean excitation potential of the atoms of the 
stopping material; it takes account of the effect of the electron binding 
energies on the rate of energy loss. It can be estimated theoretically but 
the results are not accurate enough to be useful, and in practice it must 
be determined by experiment. (See below.) The quantity 6/2 inside 
the square brackets has been added to the original Bethe formula to 
take account of the density effect discussed below. 

At very low energies the moving particle can pick up and lose electrons 
in its passage through matter, with the result that the average charge 
on the particle is reduced. At these energies Eq. 2.1 no longer 
applies. 

When the moving particle is an electron the stopping-power equation 
is modified in two ways: (1) Since in a collision with an atomic electron 
the masses of the two particles are the same, the reduced mass of the 
system of colliding particles is less and the average energy transferred 
in the collision is less. (2) Since the two particles involved in the col- 
lision are identical, for all practical purposes the one emerging from the 
collision with the greater energy must be considered as the primary 
particle, and the lower-energy one as a “knock-on” electron. Thus 
energy losses of greater than half the energy of the particle are ruled out. 
The energy-loss formula for electrons, as derived by Bethe (2), is 


dE 2re* Z mE — 
~ = Ne 5 |b ae (2 1— ps? — 1 + B*) In2 


+18 +2 — V1 — g?)? — 5 | erg/om (2.2) 


where E is the kinetic energy of the electron. At high energies the rate 
of energy loss of electrons is about 10% less than for protons of the same 
velocity. 
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Since the specific energy loss is proportional to the density, it is often 
convenient to work with the mass stopping power, defined by the equa- 


tion 
mS = p_'(—dE/dz) (2.3) 


which is independent of density. 
Another quantity that is sometimes useful is the stopping power per 
electron, which is equal to the specific energy loss divided by the electron 


density: _dB/ax 
~ Nop(Z/A) 


The stopping power per electron does not depend explicitly on the 
atomic number, Z. The only dependence on the nature of the stopping 
material is through the mean ionization potential, 7, which appears in 
the logarithmic term. Because of this logarithmic dependence on J, the 
stopping power per electron is only a slowly varying function of the 
atomic number. 

Determination of I. In Eqs. 2.1 and 2.2 all the parameters are 
well-established constants except J, the mean excitation potential of the 
atoms of the stopping material. This quantity must be determined from 
experiment before the stopping-power equations can be used. 

Bloch (4) showed that on the basis of the Fermi-Thomas model of 
the atom J should be proportional to the atomic number Z. One should 
thus be able to write J = kZ, where k is a constant that can be deter- 
mined empirically. The experimental values of k lie between 8 and 14. 

An experimental J value is obtained by measuring the energy lost 
by a particle in a known thickness of the material of interest. The loss 
can be computed by integrating the stopping-power formula between 
the appropriate limits, and a value of J is chosen to make the calculated 
loss equal to the measured one. Because electrons are so strongly scat- 
tered in matter, the measurements must be carried out with heavy 
particles. Thus Bethe (5) determined a value of 80.5 ev for air from the 
known energies and the difference between the ranges of Th C’ long- 
range a-particles and Po a-particles. Values of J for a number of solid 
materials have been determined with high-energy protons (6, 7, 8, 9, 10). 
Although there are some discrepancies in the data of different authors, 
the results of Bakker and Segre (7) have been widely accepted, and ap- 
pear to be consistent with measurements involving secondary electrons 
set in motion by y-rays (11), at least for light elements. The values ob- 
tained by Bakker and Segre are listed in Table 2.1, along with the cor- 
responding values of k. (The value for air is that found by Bethe.) It 
is evident that the “‘constant’’ k varies considerably from element to 


(2.4) 
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element. The general features of this variation have been explained 
theoretically by Brandt (10). 


Table 2.1 
EXPERIMENTAL VALUES OF THE Mean Excitation Porentiau, I 
Material I k=I1/Z 
H 15.6 15.6 
Li 34.0 11.3 
Be 60.4 15.1 
C 76.4 12.7 
Al 150 11.5 
Fe 243 9.3 
Cu 279 9.6 
Ag 428 9.1 
Pb 758 9.2 

Air 80.5 


Because the stopping power varies only as the logarithm of J, the dis- 
crepancies in the different measured J values result in much smaller dis- 
crepancies in the stopping powers calculated from them. With the 
information at present available the latter can be computed to better 
than 1% over a wide range of energies. 

Nonparticipating electrons. At low velocities, where the maximum 
energy that can be transferred to an electron (2mv? for a heavy particle 
moving with velocity v) becomes comparable to the binding energy of 
the K-electrons, the K-electrons do not contribute fully to the stopping 
power, and a correction term for these ‘“‘nonparticipating electrons”’ 
must be subtracted from the expression for the stopping power. This 
correction is discussed by Bethe and Ashkin (3), who also provide the 
data necessary for its calculation. For electrons, which have much 
higher velocities than heavy particles of the same energy, the correction 
rarely has to be considered. 

The density effect (polarization effect). Bethe’s original stopping- 
power equations were based on the assumption that the interaction of 
the moving particle with any particular atom was not affected by the 
presence of other atoms. It was later pointed out (12), however, that 
interactions with distant atoms should be influenced by intervening 
atoms: the latter will polarize in the electric field of the moving particle 
and reduce the field at a distance, with the result that energy losses by 
distant collisions will be reduced. Since the Lorentz contraction of the 
electric field of the moving particle makes distant collisions increasingly 
important at higher energies, the effect of polarization is more important 
at high energies. The magnitude of the effect will clearly depend also 
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on the density of the medium—.e., on the number of atoms polarizing 
per cm*—and for this reason it is frequently referred to as the density 
effect. It is important in solids and liquids but, except at very high 
energies, negligible in gases. 

The most complete calculations on the density effect appear in a series 
of papers by Sternheimer (13), which cover a wide range of energies and 
materials. The effect can conveniently be divided into two parts: 

(1) The zero-energy density effect, which is present even at very low 
energies and which does not depend on energy. This effect is auto- 
matically included in values of J determined from measurements in 
solids; these J values presumably differ somewhat from those that would 
be determined from measurements in the corresponding gases. As an 
example, the zero-energy effect amounts to between 1 and 2% in graph- 
ite. 

(2) The finite-energy density effect, which sets in at energies of the 
order of the rest energy of the moving particle, and increases with in- 
creasing energy. The magnitude of this effect, for the case of electrons, 
is shown in Fig. 2.1. 





0 10 20 30 
E, Mev 


Fig. 2.1. Per cent reduction in stopping power due to finite-energy density effect 
(for electrons). 
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In Eqs. 2.1 and 2.2 the symbol 6 represents the finite-energy density 
effect as calculated in the papers by Sternheimer. The zero-energy ef- 
fect in solids is generally taken into account by using J values determined 
in solids. 

Variation of stopping power with energy. The way in which the 
stopping power varies with energy can be seen from Fig. 2.2. The 
descending portions of the curves as the energy increases correspond to 
the variation as 1/v? predicted by classical physics. At energies ap- 
proaching the rest energy of the moving particle relativistic effects 
cause the stopping power to level off to a broad minimum and then 
increase gradually. In the case of liquids and solids this increase is 
modified downward by the density effect; this is illustrated by the curves 
for graphite and air, which are almost identical below 1 Mev but separate 
above that energy, the curve for the solid (graphite) falling below that 
for the gas (air). 

Linear energy transfer. A quantity closely allied to stopping power 
is the linear energy transfer, abbreviated LET. This is defined (14) as 
follows: 


Innear energy transfer (LET) is the linear-rate of loss of energy 
(locally absorbed) by an ionizing particle traversing a material 
medium. 


Linear energy transfer may be conveniently expressed in kilo elec- 
tron volts per micron (kev/y). 


The use of the concept of LET concentrates attention on the transfer 
of energy to the material, rather than on the loss by the particle. This 
is significant when the biological effect of the radiation is being con- 
sidered, since biological effect. has been found to depend on the way in 
which the energy loss is distributed along the particle tracks (i.e., on the 
LET) as well as on the total energy absorbed per gram. 

For a charged particle of a specified energy the LET value is cal- 
culated directly from the stopping-power formula. If it is considered 
that 6-rays of more than a certain energy are not “locally absorbed,’’ 
energy transfers of more than this amount must be excluded from the 
calculation (15). Under practical conditions of irradiation, such as occur 
in a beam of x-rays, a small region in the medium is crossed by secondary 
particles with energies ranging from some maximum value down to zero. 
There is thus a distribution of LET values present, and one must either 
calculate some sort of average value (16, 17, 18, 15), or derive the com- 
plete distribution (18, 19, 15). Some typical LET distributions are 
presented in reference 14. 
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Energy loss by radiation. When a swiftly moving charged particle 
of mass m and charge z passes close to a nucleus of charge Z, it ex- 
periences & Coulomb force proportional to zZ and consequently under- 
goes an acceleration proportional to 2Z/m. According to classical 
physics, an accelerated charge radiates electromagnetic energy at a rate 
proportional to the square of the acceleration, i.e., in this case propor- 
tional to z?Z?/m?. The proportionality to the square of the atomic 


5 


I (hv), arbitrary units 





9) 0.5 1.0 
hv/Eo 


Fig. 2.3. Energy distribution of bremsstrahlung for various initial electron ener- 
gies, Eo. 


number means that loss of energy by radiation should be much more 
important in heavy elements than in light, whereas the inverse pro- 
portionality to the square of the mass means that light particles should 
radiate much more readily than heavy ones. In fact electrons radiate 
more than a million times more strongly than protons, and it is only with 
the radiation of electrons that one need be concerned. 

When the radiation process is treated quantum mechanically, the 
concept of the continuous emission of weak radiation is replaced by a 
small probability that the electron will emit a photon with any energy 
between zero and the full kinetic energy of the electron. When interest 
is centered on the loss of energy by the electrons, the emitted radiation 
is usually referred to as bremsstrahlung; when the radiation itself is of 
primary interest, it is often referred to as the continuous x-ray spectrum. 
Conventional x-rays are bremsstrahlung produced deliberately under 
controlled conditions. 
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The bremsstrahlung energy spectrum extends from zero to the full 
initial energy of the electrons. Typical distributions for a number of 
electron energies are shown in Fig. 2.3. 

At low energies the emitted photon can have any direction with 
respect to the incoming electron; at high energies the radiation is di- 
rected strongly in the forward direction, the average angle with the 
direction of the electron being given approximately by 


6 = Ee (2.5) 
where £o is the kinetic energy of the electron before the emission of the 
photon. 

As pointed out above, the rate of energy loss by radiation is roughly 
proportional to the square of the atomic number of the stopping ma- 
terial. It also increases rapidly with increasing energy, being negligible 
below 100 kev and rising to values greater than the rate of loss by col- 
lision at energies between 10 and 100 Mev, depending on the material. 
In the range of energies where the screening effect of the atomic electrons 
on the nuclear charge can be neglected (mc? « E « 137mc?Z~*), the 
rate of energy loss is given by 


_ (dE\ 72 tV (LY 
(=) were ett aaah 


2(E + mc?) =| 


where N is the number of atoms per cm’, and the other symbols have the 
same meanings as in Eq. 2.1. At very high energies (EF > 137mc?Z~ %), 
where the screening by the atomic electrons is complete, the rate of loss 
is given by 


(=) N(E + me?) Z(Z+ 4) ( e y (41 183 i *) / 
—({—] = ic.) ——— | => n—, + -—] erg/cm 
de Jcs m 137 \me ze 9)" 


(2.7) 


The parameter ¢ is a function of the atomic number, and ranges between 
1.1 and 1.4. Some values of ¢ are listed in Table 2.2. 

The rate of energy loss by radiation for air, water and graphite is 
shown as a function of the electron energy in Fig. 2.2. The linear rates 
of loss have been divided by the densities, so that the curves are inde- 
pendent of the density of the stopping material. 
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Table 2.2 


VALUES OF THE PARAMETER { FOR DIFFERENT ATOMIC NUMBERS 


NX 
m~, 


8 
=8B88 


The ratio of the radiation loss to the collision loss at the energy E is 
given roughly by 
(dE/dzt)rna EZ 


(dE/dz)oo1  1600mc? 


It can be seen from Eq. 2.8 that whereas in graphite (Z = 6) the radia- 
tion loss only becomes comparable to the collision loss at energies above 
100 Mev, in lead (Z = 82) the two are equal at about 10 Mev. 

The above equations give the average rate of energy loss by radiation. 
Since the loss actually occurs in a small number of interactions, each 
involving the radiation of a large amount of energy, the fluctuations 
about this average rate are great—much greater than for collision losses. 

Elastic scattering. In addition to interacting with atomic electrons, 
a moving charged particle can also interact with the Coulomb fields of 
nearby nuclei. The result is a change in the direction of motion of the 
moving particle, and the process is called elastic scattering. When the 
moving particle is a heavy one, it must pass very close to a nucleus in 
order to suffer an appreciable deflection, and consequently such deflec- 
tions are rare and scattering is relatively unimportant for such particles. 
But an electron, because of its small mass, is deflected in much more 
distant collisions. As a result elastic scattering occurs more frequently 
than energy loss by ionization or excitation, and the path of an electron 
is marked by frequent changes of direction. 

The scattering of electrons has been discussed in some detail by Bethe 
and Ashkin (3). The cross section per atom for scattering between the 
angles @ and 6 + d@ is given by 


tow eg net oe eile 0 ) as 
= _ — — + —(sin- — , 
°  Omiy? sin* 6/2 = 2 137 9 


where v is the velocity of the moving electron and the other symbols 
have their usual meanings. It is evident from this equation that the 


(2.8) 
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scattering is greater at low velocities, and that small-angle scatterings 
are much more probable than large. 

When electrons pass through any appreciable thickness of matter, 
their net deflections are the result of a large number of individual scat- 
tering processes. These net deflections are distributed according to the 
Gaussian law about the original direction. The mean square deviation 
as a result of multiple scattering in a thickness d cm is given approximately 


by 

a 4 2 

ge? = eee eee, | 4e2Nd (—) (2.10) 
p*v* my 


where N is the number of atoms per cm’, 7 is the electron momentum, 
and h is Planck’s constant divided by 2x. The thickness d must be 
great enough so that electrons will undergo a large number of scatterings 
in passing through it. Since the mean square angle of deflection is 
roughly proportional to the thickness traversed, the root mean square 
deflection is roughly proportional to the square root of this thickness. 
An idea of the magnitude of the multiple scattering in a number of 
materials can be obtained from Table 2.3. Where no value is listed, the 





Table 2.3 
RMS ANGLE OF MULTIPLE SCATTERING FOR ELECTRONS (DEGREES) 
d (cm) 
Air Water Aluminum Lead 

E eee, ee A eee 
(Mev) 1 10 100 0.01 0.1 1 0.01 0.1 1 #£4«0.01 0.1 
0.1 18 
0.3 6.0 25 21 50 
1 2.0 8.8 33 7.3 28 18 67 91 
3 0.8 3.4 13 2.8 11 6.8 26 35 130 
10 0.3 1.1 4.2 0.9 3.6 13 2.2 8 30 12 42 


thickness of material is comparable to or greater than the electron range. 

Fermi (20) has pointed out that at extreme relativistic energies the 
root-mean-square angle of scattering is given roughly by an expression 
of the form, 


Ve = kVG/E (2.11) 


When E is expressed in Mev, d in cm, and @ in radians, k has the approx- 
imate value 0.07 for air, 2 for water, 5 for aluminum, and 24 for lead. 
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Equation 2.11 is a very-high-energy approximation; it agrees with Kq. 
2.10 to within about 20% at 10 Mev and within a factor of 2 at 1 Mev. 

Penetration of charged particles through matter. As a beam of 
charged particles passes through matter, collisions reduce the energies 
of the particles and scattering processes change their directions, and the 
eventual result of both processes is to remove particles from the beam. 
The number of particles thus tends to decrease as the thickness of ma- 
terial traversed increases. The penetration of charged particles in 
matter can be examined by interposing successively thicker sheets of 
material between a source of particles and a detector, and observing 
the number of particles as a function of absorber thickness. A graph 
of the results of such a measurement is called a ‘‘number-distance curve.”’ 
The shape of such a curve depends on the kind of particles and their 
initial energy distribution, and on the nature of the absorbing material. 

A heavy particle such as a proton or a-particle can lose only a small 
fraction of its energy in any one collision with an atomic electron, so 
that such a particle slows down gradually as a result of a large number of 


x 1 





Fig. 2.4. Typical number-distance curve for heavy particles. 


small energy losses. Furthermore, since large deflections are rare, the 
path of the particle is generally very nearly straight. Consequently a 
number of heavy particles starting out with the same energy all travel 
about the same distance in matter before stopping. The number- 
distance curve in this case has the form shown in Fig. 2.4. Of course, 
while at distances such as x, the number of particles has not changed, 
the energy of the particles has decreased. 

Because of the random nature of the individual energy losses, the dis- 
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tances travelled by different particles before stopping are not exactly the 
same, but are distributed on either side of a mean value. The mean 
distance FR travelled by charged particles in matter is called the mean 
range. In terms of the mean rate of energy loss — (dE/dx) at the energy 
E, the mean range is given by 


dE 
a f gee (2.12) 
Eo ~~ (dE/dz) 


where Ep is the initial energy. Unfortunately this equation cannot be 
used for the precise calculation of ranges simply by substituting for 
— (dE/dx) the theoretical expression given by Eq. 2.1, since the latter 
is not valid at very low energies. 

The mean range can be determined from the number~distance curve. 
Differentiating the curve with respect to x gives a curve representing the 
number of particles stopping per unit interval of distance at z (the 
dashed curve in Fig. 2.4), and the position of the centroid of this curve 
gives the mean range. Moreover, since the ranges of the individual 
particles generally lie in a Gaussian distribution about the mean, the 
mean range is equal to the most probable range, which occurs at the 
maximum of the differential curve, or the point of maximum slope on the 
number—distance curve. 

In practice it is often more accurate to measure the extrapolated 
number—distance range R, obtained by drawing a tangent to the steepest 
portion of the number-distance curve and extrapolating it to cut the 
distance axis. (See Fig. 2.4.) Once the extrapolated range is de- 
termined, the corresponding mean range can be calculated (3). 

In Fig. 2.5 the mean range for protons in air is shown as a function of 
energy. Ranges of other heavy particles can be found from the proton 
range-energy curve with the aid of a simple relationship. Since the rates 
of energy loss of two different charged particles of the same velocity are 
simply proportional to the squares of their charges (except at very low 
energies), the range of a particle of charge z and mass M is related to 
the range of a proton (charge 1, mass M,) of the same velocity by 

Rim 3 ome R 2.13 
2M “3 22 M, ?p ( : ) 
For example, a 2-Mev deuteron, which has the same velocity as a 1-Mev 
proton, has a range equal to twice that of a 1-Mev proton. Extensive 
range-energy curves for heavy particles are given by Bethe and Ashkin 
(3) and by Birge et al. (21). 

For electrons the slowing-down process is more complicated. An 

electron can lose up to half its energy in a single collision, bremsstrahlung 
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Fig. 2.5. Range—energy curve for protons in air. 


is more probable and can involve large losses, and scattering leads to a 
much more erratic path. The result is that a group of electrons starting 
out from the same point in the same direction and with the same energy 
will come to rest at widely separated points. Those that chance to travel 
in relatively straight lines and to escape large individual energy losses 
will travel a long way; those that chance to undergo many large changes 
of direction or to suffer large energy losses will come to rest nearer the 
point of origin. This situation is reflected in the number—distance curve 
for monoenergetic electrons, which has the general form shown in Fig. 
2.6. The curve falls off more or less linearly with increasing thickness, 
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and extrapolation of the linear portion to the z-axis gives the extrapolated 
range, the quantity that is commonly determined for monoenergetic 
electrons. 

When the electrons start out with a continuous distribution of energies, 
as do the §-rays from a radioactive substance, the number—distance 
curve can be thought of as arising from the superposition of a large 
number of linear curves of different slopes, corresponding to the different 
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Fig. 2.6. Typical number-distance curve for monoenergetic electrons. (1-Mev elec- 
trons in aluminum.) 


initial energies. The resulting curve is roughly exponential, as illus- 
trated by Fig. 2.7, where the logarithm of the number of particles trans- 
mitted is plotted against absorber thickness. The observed counting 
rate in the detector usually levels off to a nearly constant value at large 
thicknesses because of the presence of y-rays which are much more pene- 
trating than 6-rays. This y-ray “‘tail’’ can be extrapolated back to show 
the contribution of y-rays at smaller thicknesses. If the y-ray contribu- 
tion is subtracted from the measured curve, what is left should be due 
to the §-rays alone. 

The quantity usually determined for §-rays is the maximum range, 
which is defined as the absorber thickness at which the total curve joins 
the y-ray tail. (See Fig. 2.7.) The direct determination of the maximum 
range involves considerable uncertainty, and a number of special tech- 
niques have been devised to reduce this uncertainty. Some of these are 
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outlined in a review article by Katz and Penfold (22), in which the 
authors have also collected much of the data on electron ranges. They 
find that between 0.8 and 3.0 Mev the experimental results in aluminum 
for both the extrapolated ranges of monoenergetic electrons and the 


log N 





Fig. 2.7. Typical number-distance curve for §-rays. (1-Mev §-rays in aluminum.) 


maximum ranges of f-rays can be closely represented by an empirical 


formula: 
R = 527E — 112 (2.14) 


where R is the range in mg/cm? and E is the electron energy in Mev (or 
the end-point energy in the case of a continuous f-spectrum). Below 
0.8 Mev a somewhat more complicated formula (22) gives a closer fit to 
experiment. 

Ionization. The energy lost by a moving charged particle in col- 
lisions with the atoms of matter goes partly to excitation of these atoms 
and partly to ionization. The path of a moving charged particle is 
marked by a trail of ionized atoms. 

It is sometimes useful to distinguish two types of ionization. The 
number of ion pairs produced by the direct interactions of the moving 
particle (i.e., the number of electrons ejected by it) is called the primary 
ionization. Some of the ejected electrons—the 6-rays—have enough 
energy to produce further ionization in their turn, and this is called sec- 
ondary ionization. The two together constitute the total zonization. 
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Specific ionization refers to the number of ion pairs per unit length of 
track. 

The primary and total ionization may be examined separately in a 
cloud chamber. If the chamber is expanded just before a particle passes 
through, the droplets form on the ions and clumps of ions due to primary 
ionizing events, and a count of these droplets yields the primary specific 
ionization. On the other hand, if the chamber is not expanded until all 
the ion pairs have had a chance to diffuse apart, a droplet forms on each 
individual ion, and a count yields the total specific ionization. The 
latter is usually about three times the primary ionization. 

Figure 2.8 shows the total specific ionization of an a-particle in air as 
a function of the distance from the end of its range. (Such a plot is 
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Fig. 2.8. Bragg curves for protons and a-particles. 


known as a Bragg curve.) The specific ionization follows the rate of 
energy loss, increasing as the particle slows down, until the point is 
reached where the particle is moving so slowly that it can capture elec- 
trons and neutralize its charge. For comparison the Bragg curve for a 
proton is included in the same figure. 

If one measures the energy lost by a moving charged particle in a gas 
and the number of ion pairs produced as a result, the ratio of the two 
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quantities gives the important parameter W, the mean energy loss per 
ion pair. Since some energy goes to excitation as well as to ionization, 
W is always greater than the ionization potential of the gas. For ex- 
ample, the ionization potential of nitrogen is 15.5 ev, whereas W for 
nitrogen is approximately 35 ev per ion pair. Thus in this case about 
half the energy lost by the particle goes to ionization and half to excita- 
tion. 

A notable fact about the mean energy loss per ion pair is that it is very 
nearly independent of the energy of the moving particle, and is prac- 
tically the same for electrons as for heavy particles. Values for W in a 
number of gases are presented in Table 2.4. The value of 34.0 ev for air 


Table 2.4 
EXPERIMENTAL VALUES OF W FoR GaSsEs (24) 


(The experimental uncertainty is about 2% for the electron results and about 1% 
for the a-particle results) 


W (ev/ion pair) 


Electrons a-particles 


H2 36.3 36.5 
He 41 44 

Ne 34.7 36.5 
Os 31.0 32.4 
Air 34.0 35.3 
A 26.4 26.4 
CO: 32.8 34.4 
CH, 26.4 28 .0 
CH, 27.0 29.3 


is an average of several recent determinations (24); it coincides with the 
value recommended in reference 14. 


ELECTROMAGNETIC RADIATION 


Unlike charged particles, which tend to lose energy gradually through 
a large number of small transfers, an x-ray or y-ray tends to lose much 
or all of its energy in a single interaction. These interactions are of 
three principal kinds: the photoelectric effect, which predominates at 
low energies; the Compton effect, which predominates at intermediate 
energies (0.05 to 15 Mev in aluminum, 0.5 to 5 Mev in lead); and pair 
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production, which predominates at high energies. It will be convenient 
to discuss each of these three types of interaction separately before pro- 
ceeding to examine how they combine to produce the observed phe- 
nomena of x-ray attenuation and scattering in matter. 

The photoelectric effect. In a photoelectric interaction the entire 
energy of the photon is transferred to an atomic electron. The electron 
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Fig. 2.9. Angular distribution of photoelectrons. (Relative number per unit inter- 


val of angle.) 


emerges with a kinetic energy, E, equal to the photon energy, hv, less 
the binding energy of the electron in the atom, ¢: 


E=h-—¢ (2.15) 


The process cannot occur if the photon energy is less than the binding 
energy of the electron, nor can it occur with a completely free electron, 
since without the recoil of the rest of the atom, energy and momentum 
cannot both be conserved. Consequently the cross section for a photo- 
electric interaction with the electrons of a particular atomic shell is zero 
for photon energies less than the binding energy of an electron in that 
shell, jumps to a high value for an energy just equal to the binding 
energy, and then falls off toward higher energies. 

For photon energies greater than the K-shell binding energy, about 
80% of the photoelectric interactions are with K-electrons. In this 
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energy region the total photoelectric cross section per atom is given 
roughly by 
at = kZ*/(hv)? (2.16) 


where k is a constant and Z is the atomic number of the atoms involved. 
Equation 2.16 indicates the main features of the photoelectric cross sec- 
tion: a strong dependence on atomic number, and an almost equally 
strong dependence on photon energy. The effect is most probable at low 
energies and in high-atomic-number materials, and becomes less probable 
as the energy increases or the atomic number decreases. Exact values 
of the cross section as a function of energy in a number of materials are 
given by Bethe and Ashkin (3) and by Grodstein (25), who also present 
discussions of the theory. 

Figure 2.9 shows how the angular distnbution of the photoelectrons 
depends on the photon energy (26). At very low energies the electrons 
are ejected predominantly at right angles to the direction of the incoming 
photon, but as the energy increases the distribution shifts increasingly 
into the forward direction. 

The Compton effect. In a Compton interaction a photon of energy 
hy undergoes an elastic collision with a free, or nearly free, electron. 
The electron recoils with an energy FE at an angle ¢, and a scattered 
photon of reduced energy hy’ goes off at an angle @. (See Fig. 2.10.) 
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Fig. 2.10. The Compton effect. 


From the equations for the conservation of energy and momentum 
in the collision it can be shown that the change in wavelength in the col- 
lision is given by 


h 
’ — A = — (1 — cos 8) (2.17) 
me 


when h is Planck’s constant, m is the rest mass of the electron, and c is 
the velocity of light. Since the frequency, »v, is equal to c/X, it follows 
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that the energy of the scattered photon is 
hy 
1+ (hy/me?)(1 — cos 6) 


Equation 2.17 shows that, for a given angle of scatter, the change in 
wavelength is independent of the energy of the primary photon. Conse- 
quently the change in photon energy increases with increasing energy 
and is negligibly small at very low energies. 

The energy of a photon scattered back at 180° is 


hy 
DN ny Ne a ee 
hv’ min 1+ (hv /mey (2.19) 
and as the primary energy, hv, increases, this approaches a limiting value 
of 0.25 Mev. Since a photon scattered through a very small angle loses 
very little energy, it is evident that at primary energies much greater 
than 0.25 Mev the energy of the scattered photon depends strongly on 
the angle of scattering. 
The energy of the recoil electron is equal to the difference in energy 
between the primary and scattered photons: 


E = hv — hy’ (2.20) 


The cross section for Compton scattering is given by the Klein- 
Nishina formula (27). This formula gives for the cross section per elec- 
tron, d .c, for the scattering of a photon of energy hy into the element of 
solid angle dQ at the angle @ (with the energy hv’ as given by Eq. 2.18), 


d e (™) (- Bs a) aa (2.21) 
o = —(—) |(— + — — sin : 
4 m*c* \ hy y’ hy 


The total cross section per electron, .c, for scattering through any angle, 

is obtained by integrating Eq. 2.21 over all possible values of 6. The 

result is 

2re* - +a [= + a) 
. 1+ 2a 


a 


hy’ (2.18) 








oc >= 


1 
—-In(1+ 2a) | 


m*c* 
1+ 3a 


@ 4 2a)? + =a} (2.22) 


+ =n (1 +22) - 


where a = hy/mce’. 

The total Compton cross section .¢ is equal to the fraction of the pho- 
tons of energy hy that would undergo Compton scattering in a layer of 
material containing one electron per cm?. Of the energy removed from 
the primary beam in this way some is converted to the kinetic energy of 
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the recoil electrons and some appears as scattered photons. The fraction 
of the energy in the incoming beam that is converted to scattered 
photons, per electron/cm?, is given by the Compton scattering cross 
section, e@s, 


* 2re* 1 In (1 -+ 2a) + (1 + a)(2a? — 2a — 1) 4a? 
oe rec | 208 (1 + 2a)? 3(1 + 2a) 
(2.23) 


The fraction of the incoming energy given up to recoil electrons is then 
given by the Compton absorption cross section, .oqa, where 


eFa = eF ~~ es (2.24) 


The way in which the various Compton cross sections vary with 
photon energy is shown in Fig. 2.11. The total cross section falls off 
steadily with increasing energy, as does the scattering cross section. 
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Fig. 2.11. Variation of Compton cross sections with energy. 


On the other hand the absorption cross section increases to a maximum 
around 0.5 Mev before falling off with the other cross sections, reflecting 
the fact that the average fraction of the total photon energy going to the 
recoil electron increases with increasing energy. 
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It is sometimes useful to deal with the cross section per atom, go. This 
is just equal to the cross section per electron multiplied by the number of 
electrons per atom, i.e., 

at = Zo (2.25) 


The angular distribution of the scattered radiation is shown in Fig. 
2.12. The radial distance represents the energy scattered through an 
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Fig. 2.12. Angular distribution of Compton scattered photons. (Relative energy 
per unit interval of angle.) 


angle 6. At very low energies the scattering is symmetrical about 90°, 
and as the energy increases the scattering tends to be more and more in 
the forward direction. Figure 2.13 shows the number of recoil electrons 
ejected as a function of the angle to the direction of the incoming photon. 
There are, of course, no recoils at angles greater than 90°, and again the 
distribution becomes more strongly forward as the energy increases. 

Finally, Fig. 2.14 shows the distribution in energy of the recoil elec- 
trons for a photon energy of 1 Mev. This distribution is typical of those 
at all energies, being fairly constant at the lower energies and ending in a 
sharp rise near the maximum possible recoil energy. 

Extensive curves of angular and energy distributions for the Compton 
effect are given by Nelms (28). A number of energy distributions have 
also been tabulated by Johns et al. (29). 
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Fig. 2.14. Energy distributions of Compton recoil electrons. (hy = 1 Mev.) 
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Pair production. The third major type of x-ray interaction is pair 
production. In this process the photon disappears and is replaced by a 
pair of particles—a positron and an electron. The energy of the photon 
goes partly to the rest masses of the two particles and partly to the 
kinetic energies E, and E,. The sum of the two kinetic energies must be 
equal to the energy in excess of the total rest energy: 


E, + E, = hy — 2me? (2.26) 


Obviously the process cannot occur for photon energies less than 2mc? = 
1.02 Mev. The available kinetic energy can be divided in any propor- 
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Fig. 2.15. Energy dependence of pair-production cross section. 


tion between the two particles, and in fact all proportions are about 
equally likely. 

Pair production can only occur in the Coulomb field of a nucleus, 
since the recoil of the nucleus is necessary to conserve momentum. The 
kinetic energy of the recoiling nucleus is negligible. 

The rest energy of the two particles is removed only temporarily from 
the radiation field. The electron simply dissipates its kinetic energy and 
stops; but the positron, after losing its kinetic energy and coming almost 
to rest, eventually encounters a negative electron and the two annihilate 
with the emission, in opposite directions, of two 0.51 Mev y-rays. It is 
also possible for the positron to annihilate before losing all its kinetic 
energy, but this is much less likely. Annihilation with the emission of 
three photons can also occur with a small probability. 
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The total cross section for pair production is zero below 1.02 Mev and 
increases above that energy in the manner shown in Fig. 2.15. The cross 
section is very nearly proportional to Z?, so that the effect is much more 
important in heavy materials than in light. What is plotted in Fig. 2.15 
is the cross section divided by Z?, for the cases of carbon and lead. The 
curves for elements of intervening atomic numbers would lie between 
these two extreme curves. 

The angular distribution of the positrons and electrons follows the 
same general trend with energy as the distributions of photoelectrons 
and Compton recoils, tending more toward the forward direction as the 
energy increases. 

Photonuclear reactions. One further type of interaction between 
x-rays and matter bears mentioning. At energies above the binding 
energy of a nucleon in a nucleus (of the order of 10 Mev for most ele- 
ments) it is possible for a photon to eject a neutron or proton from a 
nucleus. In such a reaction the energy of the photon is completely 
absorbed, part going to separate the nucleon from the nucleus and the 
rest to kinetic energy. The cross section for photonuclear reactions is 
seldom more than a few per cent of the total cross section for interactions 
with atomic electrons. The cross section always has a threshold energy 
below which the reaction does not take place, and generally rises to a 
maximum with increasing energy and then falls off again. Photonuclear 
reactions are seldom important in dosimetry. 

The attenuation of x-rays in matter. In many experimental 
situations what is of concern is not the probability that an x-ray photon 
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Fig. 2.16. Narrow-beam geometry. 


will undergo an interaction of a specific type, but the attenuation of a 
beam of x-rays due to the combined effects of all possible interactions. 
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Consider a narrow beam of x-rays falling on a detector which measures 
the intensity, J. (See Fig. 2.16.) If a small thickness dz of absorbing 
material is interposed in the beam, the intensity will be reduced by some 
amount d/J, as a result of the interactions of some of the photons with 
atoms of the material. The fraction of the beam undergoing photoelec- 
tric interactions will be equal to the cross section per atom, gr, multiplied 
by the number of atoms per cm?, pNo dx/A, p being the density in g/cm®, 
No Avogadro’s number, and A the atomic weight. Thus the reduction 
in the intensity due to photoelectric interactions will be 


nf 
dl, = —gr = I dz (2.27) 
Similarly the reduction due to Compton interaction will be 
pNo 
dI, = “eer ae (2.28) 


where the total Compton cross section is used, because in a Compton 
interaction all the energy of the photon is removed from the beam—some 
in the kinetic energy of the recoil electron and some in the scattered 
photon which, under the ‘‘narrow-beam’”’ conditions envisaged here, goes 
off in a different direction and misses the detector. Finally the reduction 
due to pair production is 


N, 
dl, = —x —"I dz (2.29) 
A 
The total reduction in intensity is given by 


N 
i = dl, dha = — (ar + ae tat) Idx (2.30) 


or 
dI = —pl dz (2.31) 
where 
pNo 
ae (at + ao + ak) (2.32) 


The quantity » defined by Eq. 2.32 is called the linear absorption coef- 
ficient (or linear attenuation coefficient) ; it is a function of the absorbing 
material and photon energy but not of the radiation intensity or absorber 
thickness. It can be thought of as made up of partial coefficients for 
each of the various types of interaction: 


p=ertote« (2.33) 


pNo 
where r = —— gr, etc. 


36 PRINCIPLES OF RADIATION DOSIMETRY 


The relationship between the intensity at the detector and the thick- 
ness, x, of absorber placed in the beam is obtained by integrating Eq. 
2.31: 

I = Ine (2.34) 


where [ is the intensity with no absorber present. Equation 2.34 repre- 
sents the familiar exponential law of absorption observed when attenua- 
tion measurements are made on photons of a single energy under narrow- 
beam conditions, 1.e., conditions such that a negligible fraction of the 
photons scattered in the absorber can reach the detector. 
When Eq. 2.34 is to be represented graphically it is usually expressed 
in the form 
InI =InIo — px (2.35) 


A plot of In J vs. x gives a straight line of negative slope. 

It is sometimes useful to characterize the absorption of radiation in a 
particular material by the half-value layer (HVL), which is defined as the 
thickness of absorber that will reduce the intensity of a narrow beam of 
the radiation to one half. When the absorption is exponential, the half- 
value layer is 

HVL = In 2/nu = 0.693/u (2.36) 


The simple exponential absorption described by Eq. 2.34 is observed 
only for a mono-energetic beam of x-rays. When photons of more than 
one energy are present, the net absorption coefficient is determined by 
the spectrum of photon energies. As the absorber thickness is increased 
under narrow-beam conditions, the components with the highest absorp- 
tion coefficients are removed faster than the others and become progres- 
sively less important. Consequently the effective absorption coefficient 
becomes smaller as the absorber thickness is increased; one speaks of 
the beam becoming “harder.’”’? On a semilogarithmic plot of the type 
mentioned above, the attenuation curve, instead of being a straight 
line, shows a continuously decreasing slope. 

The linear absorption coefficient, u, is proportional to the density of 
the absorbing material. It is often more convenient to deal with the mass 


absorption coefficient, u/p, 
N 
n/p = oe Gps thaw (2.37) 


which is independent of density. In terms of the mass absorption coef- 


ficient, Eq. 2.34 becomes 
I = Ige  loez (2.38) 


where pz is the absorber thickness in g/cm?. 
The mass absorption coefficient is particularly useful in the range of 
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energies (and atomic numbers) where the Compton effect is the only 
significant interaction, for then 
No N Z 
n/p = o/p = A ao = ar al (2.39) 
Since .o does not depend on atomic number and Z/A is close to 0.5 for 
all low-atomic-number materials, u/p is very nearly the same for all 
light materials over this energy range. 

Extensive tables of x-ray absorption coefficients have been calculated 
by Davisson and Evans (26) and more recently by Grodstein (25). In 
both publications the agreement between the theoretical coefficients and 
experiment is discussed. Over most of the range of energies and atomic 
numbers examined, the agreement is to better than 1%. 

It was mentioned above that in an experimental arrangement of the 
sort shown in Fig. 2.16 only the radiation emerging from the absorber 
in the original direction of the beam is detected, and what is measured 
is the total absorption coefficient. However, not all of the energy that 
fails to reach the detector when an absorber is placed in the beam is 
actually retained in the absorber; some escapes in other directions in the 
form of scattered photons. The coefficient yu, that would be measured 
if all this scattered radiation were detected, i.e., the coefficient referring 
to the actual absorption of energy in the absorber, is called the energy 
absorption coefficient or true absorption coefficient. Like the total coef- 
ficient, the energy coefficient is the sum of three partial coefficients: 





Ma = Ta + Oa + Ka (2.40) 
Here 7, 1s the coefficient for photoelectric absorption, 
pce 
T= r( . *) (2.41) 
hy 


where f is the probability that the binding energy, ¢, of the electron 
ejected from the atom will be re-emitted as a characteristic x-ray; at 
photon energies large compared to the K-shell binding energy the factor 
in parentheses is very nearly equal to unity and 7, is, for all practical 
purposes, equal to r. The quantity o, is the Compton absorption co- 
efficient, given by 


| Z 
Ca = pNo A eF%a (2.42) 
Finally, the pair-production coefficient is 
hy — 2mc? 
Kg = x | ————_ (2.43) 
hy 


where the factor in parentheses is the fraction of the photon energy which 
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appears as kinetic energy of the positron-electron pair and hence is 
absorbed near the point of interaction. The energy 2mc* reappears as 
two 0.51 Mev photons when the positron annihilates. 
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Fig. 2.17. Mass absorption coefficients for air. 
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The quantity yuo/p is called the energy, or true, mass absorption co- 
efficient. Values of this coefficient for a number of energies and materials 
are tabulated in the 1956 report of the ICRU (14) and in a paper by 
Fano (30). 

The relationship between the various absorption coefficients discussed 
above is shown in Fig. 2.17, for the case of air. Except for the extreme 
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low-energy and high-energy ends of the range, the curves for other light 
substances such as carbon, water and aluminum do not differ from those 
for air by more than a few per cent. The figure shows clearly the fact 
that between about 0.1 Mev and 2 Mev the only significant interaction 
in air is the Compton effect. Within this range of energies, sometimes 
referred to as the Compton region, the total absorption coefficient is 
equal to o and the energy absorption coefficient is equal to ag. 
Broad-beam attenuation. The reduction in the response of a de- 
tector that will be produced by any thickness of absorber in an arrange- 
ment like that of Fig. 2.16 can be predicted directly from the theoretical 
total absorption coefficient with the aid of Eq. 2.34. However, if the 
experimental conditions are less restricted (so that some of the radiation 
scattered in the absorber can reach the detector) the situation is con- 
siderably more complicated. For example, in the arrangement shown in 


Absorber 





Fig. 2.18. Broad-beam geometry. 


Fig. 2.18, the loss of photons such as A by scattering away from the de- 
tector will be partially balanced by the scattering of photons such as B 
toward the detector. Thus the effect of the absorber is to reduce the 
reading of the detector by less than it would under narrow-beam condi- 
tions. In these circumstances one can define an “effective” absorption 
coefficient to describe the attenuation over any limited range of thick- 
nesses, but this coefficient will actually be a function of the thickness. 
A more common procedure has been to replace Eq. 2.34 by one of the 
form 
T= I)Be-** (2.44) 
where the build-up factor B takes into account the addition to the ob- 
served intensity due to scattered radiation. In the most general case, 
B is equal to the ratio of the total detector response at the point of in- 
terest to the response due to the primary (i.e., unscattered) radiation 
at that point. The build-up factor depends on the nature of the detector 
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and the geometrical arrangement, as well as on the kind and thickness of 
absorber and the primary photon energy. The calculation of build-up 
factors presents a complex problem. An outline of the considerations 
involved has been given by Fano (31), who also tabulates build-up factors 
for a point isotropic source in an extended medium. More extensive 
numerical results have been published by Goldstein and Wilkins (32). 
These show that the radiation on one side of a thick barrier due to the 
presence of a source on the other side may be several times what one 
would predict on the basis of the narrow-beam equation. 

Backscattering. Another consequence of the scattering of radiation 
in matter is that when a beam of radiation strikes the surface of a solid 
or liquid medium, the spherical intensity in front of the surface is greater 
than it would have been in the absence of the medium. The additional 
intensity arises from radiation scattered through angles of more than 
90° in the medium and returning through the surface. The amount of 
this backscattered radiation depends on the photon energy of the incident 
radiation, the area of the beam, and the atomic number of the medium. 
Depending on all these factors the increase may be as much as 50%, or 
as little as a fraction of a per cent. 


NEUTRONS 


Neutrons, since they are electrically neutral, do not lose energy by 
ionization. They interact almost exclusively with atomic nuclei, the 
most common types of interaction being elastic scattering, inelastic 
scattering, nuclear reactions, and capture (33). 

Elastic scattering. The most probable interaction for a neutron of 
several Mev energy is elastic scattering. This may be thought of as 
a simple “billiard-ball” collision, in which the energy transferred and 
the angle of deflection are related by elementary mechanics. The maxi- 
mum energy that can be transferred to a nucleus of atomic weight A 
by a neutron of energy Ep (and atomic weight 1) is given by 


(=) _ A pe 
Eo/ max (A +1)? ie 


This is greatest for hydrogen (A = 1) and decreases with increasing 
atomic weight. Below about 10 Mev in hydrogen, and at lower energies 
in other elements, the scattering is isotropic in the center-of-mass system. 
It follows that energy losses of all sizes between zero and the maximum 
are equally probable. The energy lost by the neutron appears, of course, 
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as kinetic energy of the recoiling nucleus, and is quickly dissipated by 
ionization and excitation of nearby atoms. 

A useful parameter in connection with the slowing down of neutrons 
by elastic collisions is the average energy decrement, £, defined as the 
average change in the natural logarithm of the energy per collision. 
In terms of the atomic weight of the nuclei in the slowing-down material, 


£ is given by 
A — 1)? A+1 
Pa <5 n() (2.46) 
2A — 1 


The larger the atomic number, the smaller is the average energy decre- 
ment. (The formula breaks down in the case of hydrogen, for which, 
however, it can be shown that — = 1.) The average number of collisions 
to reduce a neutron from an initial energy Eo to a lower energy E is given 
by 

_ _ B Go/B) 


g 


Thus the average number of collisions to reduce the energy of a neutron 
from 2 Mev to 0.025 ev (thermal energy) in hydrogen is 18, in graphite 
over 100, and in lead about 2000. 

Scattering cross sections tend to increase with decreasing neutron 
energy, and hence a neutron slowing down in matter by elastic collisions 
tends to travel progressively shorter distances between collisions as its 
energy is reduced. 

Inelastic scattering. It sometimes happens that when a neutron 
collides with a nucleus the nucleus is left with some additional excitation 
energy and the energy of the neutron is correspondingly reduced. Such 
encounters are known as inelastic scattering. Unlike the recoil energy of 
the nucleus, the excitation energy does not appear nearby as ionization; 
instead it is eventually emitted as gamma radiation which may be ab- 
sorbed some distance away or escape from the medium altogether. At 
energies below several Mev the cross section for inelastic scattering is 
small compared to that for elastic scattering, and at energies below that 
of the lowest excited state of the nucleus (generally a few hundred kev) 
the process cannot take place at all. 

Nuclear reactions. The third possibility when a neutron collides 
with a nucleus is that it will undergo a nuclear reaction. In this case 
the neutron is incorporated into the nucleus and another particle such as 
a proton or a-particle is emitted. The kinetic energy of the reaction prod- 
ucts can be either more or less than the energy of the incoming neutron, 
depending upon the particular reaction. Except in the case of the (n, 2n) 


(2.47) 
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reaction, the emitted particles are charged and deposit their energy 
through ionization and excitation near the point of interaction. A very 
few reactions, such as the (n, a) reaction with B?°, can take place at low 
or thermal energies, but in most cases the cross section does not become 
significant below several Mev. 

Capture. At thermal energies the most probable interaction is cap- 
ture, in which the neutron is simply absorbed by the nucleus and, in most 
cases, one or more y-rays are emitted. The cross section for capture 
tends to be low at high energies, to pass through a series of resonances 
between 10 Mev and 1 ev, and then at still lower energies to increase in 
inverse proportion to the neutron velocity. 

Slowing down and capture in water. The relative importance of 
the various neutron interactions in a representative practical situation 
is illustrated by the behavior of 10-Mev neutrons in water. At all 
energies the most probable interaction is elastic scattering, and most of 
the neutrons lose energy in a series of elastic collisions until they are in 
thermal equilibrium with the nuclei of the water. This slowing-down 
process takes place in an average distance of roughly 10 cm. Since the 
scattering cross sections of hydrogen and oxygen are of the same order of 
magnitude, there are comparable numbers of collisions with nuclei of 
both elements, but, since the average energy loss in a collision with a 
hydrogen nucleus is several times greater than in a collision with an 
oxygen nucleus, most of the energy is transferred to recoil protons. 
There is no inelastic scattering in hydrogen, and, since the energy levels 
in the predominant isotope of oxygen are high and widely spaced, there 
is little inelastic scattering in oxygen except at the highest energies. A 
few per cent of the neutrons react with oxygen at energies greater than 
1 Mev; at lower energies reactions in oxygen are negligible. Once having 
reached thermal energies, the neutrons diffuse until they are captured 
by the hydrogen to form deuterium. This process is accompanied by the 
emission of a 2.2 Mev y-ray. 

Of the original kinetic energy of the 10 Mev neutrons, a small fraction 
thus goes to nuclear reactions and inelastic scattering in the oxygen and 
the rest to the recoil energy of hydrogen and oxygen nuclei. The recoil 
energy is promptly dissipated as ionization and excitation of the atoms 
of the medium. An additional 2.2 Mev per neutron is released as 
y-radiation on capture in the hydrogen. 


C. 
R. 
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CHAPTER 


D 


Determination 
of x-ray spectral 
distributions 


X-ray spectra. X-rays are generated by accelerating electrons to an 
appropriate energy and then allowing them to stop in a target composed 
of a high-atomic-number material. As discussed in Chapter 2, a fraction 
of the electron energy appears as electromagnetic radiation (brems- 
strahlung) with photon energies distributed between zero and the full 
energy of the electrons striking the target. 

The probability of the production of a photon of a particular energy 
hy decreases as the photon energy approaches the limiting value, hyo, so 
that the spectrum falls off smoothly toward zero at the high-energy end. 
In a practical x-ray machine the low-energy radiation is partially ab- 
sorbed in the target itself and in the window of the x-ray tube, and this 
absorption increases with decreasing energy. Consequently the spec- 
trum also falls off toward zero at very low energies. A typical continuous 
x-ray spectrum thus rises to a maximum at a photon energy well below 
that of the electrons producing it and then falls off to zero again, as 
shown in Fig. 3.1. (The two peaks shown in the figure are discussed 
below.) Increasing the electron energy shifts both the upper limit and 
the position of the maximum toward higher energies. Inserting an 
absorbing filter in the beam reduces the intensity at all photon energies, 
but not by the same amount: generally the reduction is greater at lower 
energies and the mean energy of the filtered spectrum is greater than 
that of the unfiltered one. 

Accompanying the continuous distribution of energy from the target 
are photons of a few discrete energies, arising from transitions of excited 
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atoms to lower energy levels. The most important of these characteristic 
lines are the K- and L-lines, which result from electrons dropping to fill 
vacancies in the K- and L-shells of atoms in the target. In tungsten, 
which is a common target material, the K-lines occur at around 70 kev 
and the L-lines at around 9 kev. The characteristic radiation shows up 
in a measured spectrum as a series of narrow peaks superimposed on the 
continuous distribution, as illustrated in Fig. 3.1. 


Intensity 
a 


hy 
Fig. 3.1. Typical x-ray spectral distribution. 


Information about x-ray spectra has been obtained in four principal 
ways: calculation, analysis of attenuation curves, crystal spectrometry, 
and scintillation spectrometry. These methods differ considerably in ac- 
curacy and in the amount of detail they can reveal. 

Calculation. In 1923 Kramers (1) published a semi-classical theory 
of the continuous x-ray spectrum in which he found that the intensity 
distribution of the radiation from a thick target, averaged over all 
directions, could be represented approximately by 


I(hv) = constant X (hyo — hv) (3.1) 


where [(hv) d(hv) is the radiation intensity in the energy range d(hv) 
at hy, and hyo is the upper limit of the photon energy distribution, de- 
termined by the tube potential. This formula has found experimental 
support at potentials below 50 kv in the work of Kulenkampff (2, 3) with 
a crystal spectrometer, although his most recent results (4), which in- 
clude a correction for the wavelength response of the spectrometer, show 
deviations from the theory. 

However well Kramers’ formula may represent the radiation leaving 
the target, the spectrum in which one is interested in practice is of course 
very different, since in passing through the tube window (and any addi- 
tional filtration) the shape of the spectrum is considerably modified. 
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A number of calculations of such filtered spectra have been carried out 
(5, 6, 7); where it has been possible to compare the results with experi- 
ment, the agreement has tended to be only approximate. (See, for 
example, reference 8.) 

Kramers’ formula applies only to very low energies. A more exact 
theory, valid at kinetic energies below the rest energy of the electron, 
has been worked out by Sommerfeld (9), and numerical values of spectral 
distributions and efficiencies of production have been calculated by Kirk- 
patrick and Wiedmann (10). There are also approximate theories ap- 
plicable at higher energies (11, 12, 13); in general these agree approxi- 
mately with measurements of spectral distributions. 

Analysis of attenuation curves. This method of determining spec- 
tral distributions is based on the following considerations. It was 
pointed out in Chapter 2 that if the transmission of a beam of x-rays of 
a single energy is measured through increasing thicknesses of some absorb- 
ing material under narrow-beam conditions, the intensity simply falls 
off exponentially with the thickness. A plot of the logarithm of the 
intensity vs. absorber thickness yields a straight line of negative slope, 
the value of this slope depending on the absorption coefficient of the 
material for that particular photon energy. However, if the beam con- 
tains a continuous spectrum of photon energies, and if the absorption 
coefficient of the absorbing material decreases with increasing energy 
(as it does at energies below several Mev), the low-energy radiation is 
more strongly absorbed than the high-energy radiation. Consequently, 
as the absorber thickness is increased, the average energy of the radia- 
tion emerging from it increases, the earlier layers having filtered out most 
of the lower-energy components. As a result the slope of the logarithmic 
transmission curve is no longer constant but becomes less steep at 
greater thicknesses. The exact shape of the transmission curve is de- 
termined by the radiation spectrum and the way in which the absorption 
coefficient varies with energy. Conversely, from the measured shape of 
the transmission curve (and the known variation of the absorption co- 
efficient with energy) it is possible to deduce the corresponding spectral 
distribution. 

The actual determination of the spectrum is carried out by fitting an 
analytical function to the transmission curve. The spectrum is obtained 
from the inverse Laplace transform of this function. The results are 
sensitive to the form of the function chosen to fit the data,* and in fact 
a function must be chosen which transforms into more or less the ex- 
pected spectral shape if incongruous results are to be avoided. The at- 


* The inevitable uncertainties in the measured results allow some leeway in the 
choice of such a function. 
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tenuation method also suffers from an inability to resolve sharp peaks in 
the spectrum, unless the energies and relative intensities of the peaks are 
already known. It works best for simple, moderately filtered spectra. 

A review of the attenuation method, including a tabulation of relevant 
absorption coefficients, has been carried out by Greening (14). 

Crystal spectrometry. Although the method of analyzing x-ray 
spectra based on the dependence of the angle of Bragg reflection in a 
crystal on the wavelength of the radiation dates back to 1912, surpris- 
ingly little work has been done on the continuous spectrum. Most 
published data have not been corrected for the dependence on wave- 
length of the reflection efficiency of the spectrometer. An exception is 
the work of Kulenkampff (3, 4), in which an attempt has been made to 
account for all the factors influencing the measured spectral shape. 
Crystal spectrometry has much the highest resolution of the techniques 
discussed here, and the method has been widely used to determine the 
energies of the characteristic lines. 

Scintillation spectrometry. This is a comparatively new tech- 
nique, the possibilities of which are still being explored (15, 16, 17, 18, 
19). The radiation is allowed to strike a sodium iodide crystal large 
enough to absorb completely most of the photons that enter it. Each 
photon that is totally absorbed in the crystal produces a flash of light 
proportional to its energy, and the light flash is converted to a voltage 
pulse with the aid of a photomultiplier. The distribution in size of these 
pulses, when corrected for the absorption of the radiation between the 
source and the crystal, the detection efficiency of the crystal, and the 
finite resolution of the spectrometer, yields the spectral distribution of 
the radiation. 

The scintillation method is most easily applied at photon energies be- 
tween about 20 and 300 kev. Below 20 kev the absorption correction 
becomes large and the resolution of the system becomes progressively 
poorer. Above about 300 kev the size of the crystal necessary to absorb 
a large enough fraction of the radiation falling on it increases rapidly. 
In a crystal of inadequate size there is an appreciable chance of a Comp- 
ton scattering followed by the escape from the crystal of the scattered 
photon; this process results in a pulse corresponding to only part of the 
energy of the incident photon, and thus tends to distort the observed 
spectrum. Increasing the size of the crystal increases the probability 
that the scattered photon will be absorbed before reaching the edge of 
the crystal. 

In many cases a knowledge of a particular spectrum can be obtained 
without actual measurement. For a number of common tube voltages 
and filtrations x-ray spectra have been measured and published; to the 
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extent that one can duplicate the original conditions of production, one 
can also expect to duplicate the measured spectrum. A number of such 
spectra are presented in reference 8, along with a list of further references 
to work at low energies. Additional data have since become available for 
tube voltages of 50 kvp (20), 30 to 100 kvp (21), and 280 kvp (22). 
Further references to work at higher energies are given by Starfelt and 
Koch (19) and Robson and Gregg (23). 

Gamma-ray spectra. Gamma-rays arise from transitions between 
definite nuclear energy levels, and consequently the y-ray spectrum of a 
radioactive substance always consists of one or more discrete lines. The 
energies and relative intensities of the y-rays from many hundreds of 
radioactive nuclides have been carefully measured and tabulated (24), 
so that in most cases to know the nuclide is to know the spectrum. How- 
ever, account must sometimes be taken of such complicating factors 
as (a) the presence in the source of traces of radioactive impurities, 
(b) bremsstrahlung resulting from the stopping of high energy §-rays 
emitted by the source, and (c) radiation degraded in energy by Compton 
scattering either in the source material itself or in the material enclosing 
it; the amount of this degraded radiation increases with the mass of ma- 
terial in and around the source. Since in the region of photon energies 
where most +-rays occur the interaction properties of the radiation vary 
only slowly with energy, the presence of a few per cent of extraneous 
radiation in the spectrum is often of minor importance in dosimetric 
problems, provided only that it 1s included in the total intensity. 
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CHAPTER 


4, 


Measurement 
of x-ray intensity 


The calorimetric method. The most direct way to measure radia- 
tion intensity is to absorb the radiation in a block of metal which is 
thermally insulated from its surroundings and measure the rate at which 
the temperature of the block increases. If all the energy striking the 
block can be converted to heat,* the intensity can be obtained from a 
knowledge of the thermal properties of the block. Although calorimetric 
techniques are difficult, principally because of the very small amounts of 
energy involved, the development of thermistors (small temperature- 
sensitive elements composed of semiconductors with high temperature 
coefficients of resistance) and the increasing emphasis on absolute meas- 
urements of energy in dosimetry have led to a recent growth in the use 
of such techniques. 

Figure 4.1 depicts the essential features of a calorimeter for measuring 
radiation intensity. A portion of the x-ray beam, defined by an aperture 
of accurately known area, passes through a thin window of aluminum or 
beryllium and strikes an absorbing block of a heavy metal such as lead 
or tungsten, in which almost all the radiation energy is turned to heat. 
The block is suspended by fine threads in an evacuated enclosure, so 
that losses of heat to the surroundings are kept small, and the enclosure 
is surrounded by a bath of water or oil kept at a constant temperature, 
so that what losses there are do not fluctuate rapidly. A thermistor 
connected to a sensitive bridge circuit detects temperature changes in 
the block, and a small electrical heating coil embedded in the block al- 
lows energy to be fed in from an external circuit. 


* The energy retained in an absorbing material does not necessarily all appear as 
heat: in some materials, for example, a fraction of it may go to the breaking of 
chemical bonds. | 
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The system is first calibrated by supplying energy to the block at a 
known rate through the heating coil and noting the response of the de- 
tector. After the system has had a chance to approach thermal equilib- 
rium again, the x-ray beam is allowed to strike the absorbing block, and 
the response of the detecting system is observed once more. From these 
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Fig. 4.1. Principal features of a calorimeter for the measurement of radiation 
intensity. 


two readings the rate at which x-ray energy was crossing the area of the 
aperture can be determined. This gives the average plane intensity 
over the area of the aperture. 

A number of factors must be taken into account in interpreting the 
results of such a measurement. At low energies (below about 50 kev) 
the attenuation of the beam in the air and in the window of the calorime- 
ter may be appreciable. Also, at the lower energies a significant fraction 
of the radiation may be backscattered from the front face of the absorb- 
ing block and hence escape detection. A correction for this effect is made 
difficult by the fact that it depends on the energy spectrum of the radia- 
tion, which is seldom known accurately. It is also necessary to take ac- 
count of any radiation that penetrates through the absorbing block or is 
scattered and escapes through the sides; losses of this kind can be esti- 
mated fairly accurately. 

At photon energies above a few tens of kev the corrections listed above 
do not present a serious problem. For example, Genna and Laughlin 
(1) have made calorimetric measurements on a beam of Co® y-rays for 
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which they quote an uncertainty of less than 2%. Other measurements 
have ranged from 400 kv to 37 Mev (2, 3, 4, 5). 

A measurement of the sort outlined above is suited to a parallel, or 
nearly parallel, beam of radiation. When the radiation is distributed 
over a wide range of directions, such a measurement will yield an approxi- 
mation to the plane intensity in a particular direction. To measure the 
spherical intensity one would require an arrangement whereby a spherical 
absorbing block could receive radiation from all directions. The techni- 
cal difficulties involved in such an arrangement have so far proved 
insurmountable. 

Calorimetric measurements are more elaborate and involved than cer- 
tain types of ionization measurements that will be described further on, 
but they have the advantage that they measure energy directly, and for 
this reason they provide a valuable check on less direct methods. 

The ionization method. Another approach to the problem of 
measuring intensity is to use the ionization produced by the radiation 
in a layer of air. Suppose, for example, that a beam of x-rays is allowed 
to pass through an air-filled gap between a pair of parallel conducting 
plates maintained at different potentials (see Fig. 4.2). The secondary 
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Fig. 4.2. Schematic diagram of a device for measuring relative x-ray intensity. 


electrons set in motion by the radiation will ionize the air between the 
plates, and the motion of the ions in the electric field will give rise to a 
current in the external circuit. Within certain limits this ionization cur- 
rent will be proportional to the intensity of the radiation, provided the 
spectrum is not changed. Such an arrangement thus yields a measure of 
relative intensity for a fixed spectrum. 

In order to make the measurement absolute, one might confine one’s 
attention to the ionization due to the interaction of the radiation in a 
fixed mass of air—say that occupying 1 cm® at 0° C and 76 cm of mer- 
cury—and attempt to measure the total ionization due to all the second- 
ary electrons originating in that mass. (This situation is illustrated in 
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Fig. 4.3.) Since the average energy required to produce an ion pair in 
air appears to be constant over a large range of electron energies (see 
Chapter 2), this ionization will be closely proportional to the energy lost 
by the radiation in the specified mass, 1.e. 


J’, = E'4/Wa (4.1) 


where J’, is the ionization produced by the secondary electrons originat- 
ing in unit mass of air, £’, is the energy lost by the radiation in unit mass 
of air, and W, is the mean energy expended by an electron per ion pair 





Fig. 4.3. X-rays ejecting electrons from 1 cm? of air. The total ionization produced 
by all the ejected electrons is related to the radiation intensity. 


produced in air. The primes are used to indicate that the corresponding 
quantities are associated with the energy lost to secondary electrons by 
the radiation at the point of interest. Unprimed quantities will be used 
in connection with the energy actually imparted to matter by secondary 
electrons at a point. 

When radiation of only one photon energy is present, the energy per 
gram E£’, absorbed from the radiation field is related to the spherical 
intensity, J,, by the expression 


E's = (uo/p) alt (4.2) 


where (u/p)a4 is the energy mass absorption coefficient for the radiation 
in air, as defined in Chapter 2, and ¢ is the time during which the radia- 
tion is absorbed. 

It follows from Eqs. 4.1 and 4.2 that 


_ ‘AWa 
* — t(ta/p) a 


This equation holds when the quantities involved are expressed in a 
consistent set of units. If, as is commonly the case, J, is expressed in 


(4.3) 
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ergs/cm?-sec, J’4 in esu/g, Wa in ev/ion pair, ¢ in sec, and (¢/p)4 in 
cm?/g, the right-hand side must be multiplied by 
1.602 X 107}? erg/ev 
4.803 X 10~?° esu/ion pair 
and Eq. 4.3 must be written 
J'aWa 


I, = 3.336 X 10~* ———— erg/cm?-sec (4.4) 
t(ua/ p)A 
It is apparent that the intensity is related to the rate of ionization 
through the absorption coefficient (ue/p)4, which is a function of the 
photon energy. (The variation of (ue/p)4 with energy is shown in Fig. 
2.17.) Thus the same rate of ionization would correspond to different 
intensities at different photon energies. 


In the case of a spectrum of energies, J, = f I,(hv) d(hv) the ioniza- 


tion is given by 


tf (uo/e)aTa(hw) a0) 
Dd 4.5 
4=—3336x103w, ve > 
whence 
I,(hv) d(h 
a 2.000) a) 
I, = 3.336 X 10-3 44 | _~ | erg/em?-sec (4.6) 


| fGwa/oale(hw) dC 


The ratio of the two integrals simply gives an average of (u,/p)4 over the 
radiation spectrum. For the purpose of evaluating this average it is 
necessary to know the relative, though not the absolute, distribution 
of intensity as a function of energy. 

How the quantity J’4 can be determined experimentally will be 
discussed in the next chapter. 


REFERENCES 


1. 8. Genna and J. 8. Laughlin, Radiology 65, 394 (1955). 

2. J. S. Laughlin, J. W. Beattie, W. J. Henderson, and R. A. Harvey, Am. J. 
Roentgenol. 70, 294 (1953). 

3. G. W. Dolphin and G. 8. Innes, Phys. Med. Biol. 1, 161 (1956). 

4. B. Zendle, H. W. Koch, J. McElhinney, and J. W. Boag, Radtatton Research 5, 
107 (1956). 

5. J. McElhinney, B. Zendle, and S. R. Domen, Radiation Research 6, 40 (1957). 


CHAPTER 


D 


Measurement 
of exposure dose 


In many instances the radiation spectrum is not known, and it is not 
possible to deduce the intensity from the observed ionization. Further- 
more, what is often of primary interest is the absorbed dose (defined in 
Chapter 1), which is more closely related to the ionization than to the 
intensity. In either case it is customary to work with a quantity called 
the exposure dose, based on the ionization in air. In 1956 exposure dose 
was defined as follows (1): 


Exposure dose of x- or gamma-radiation at a certain place is a measure 
of the radiation that is based upon its ability to produce ionization. 

The unit of exposure dose of x- or gamma-radiation is the roentgen 
(r). One roentgen is an exposure dose of x- or gamma-radiation such 
that the associated corpuscular emission per 0.001293 g of air pro- 
duces, in air, ions carrying | electrostatic unit of quantity of electricity 
of either sign. 

Exposure dose rate is the exposure dose per unit time. The unit of 
exposure dose rate is the roentgen per unit time. 


Note that 0.001293 g is just the mass of 1 cm? of dry air at 0° C and 
76 cm Hg. 

For the sake of brevity, and in order to emphasize the distinction 
between absorbed dose and exposure dose, the latter is sometimes re- 
ferred to simply as exposure.* 


* Until recently the term ‘‘dose’”’ was sometimes used to refer to what is now called 
exposure dose, and sometimes to refer to what is now called absorbed dose. Thus, 
depending on the context, the phrase ‘‘a dose of one roentgen’’ could mean either ‘‘an 
exposure dose of one roentgen’’ or “‘an absorbed dose of the magnitude that would 
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It follows from the definition of the roentgen that it is directly related 
to the quantity J’, defined in the previous chapter. In fact, if R is the 
exposure in roentgens, then 


R = 0.001293/’4 roentgens (5.1) 


where J’, is the total ionization, in esu, produced by all the secondary 
electrons set in motion by the radiation per gram of air. 

It follows from Eqs. 5.1 and 4.4 that the exposure is related to the 
quantity of radiation (intensity X time) and exposure rate is related to 
intensity. It should be noted that, like the intensity, the exposure at a 
point is simply a property of the radiation at that point; it has nothing 
to do with what happens to the radiation—whether it is absorbed or 
causes ionization in matter or passes through the region in question 
unaffected. It refers to the way in which the radiation would interact 
with air at that point, if it had the chance. 

The exposure as such is seldom of direct interest. Over a considerable 
range of energies it provides a means by which absorbed dose in tissue 
can be determined, and, in circumstances where the intensity cannot 
be measured directly, the exposure rate can provide a useful substitute. 
Because it is readily determined, the exposure is often what is obtained 
experimentally, and the other quantities are derived from it. 

The relationship between exposure and quantity of radiation or be- 
tween exposure rate and intensity, can be obtained by combining Eqs. 
5.1 and 4.4: 

R = 0.3876(ua/p) al st 
Wa 


The ratio of quantity of radiation to exposure, expressed in ergs/cm? 
per roentgen, is thus equal to W4/0.3876(ue/p)a, which depends on the 
photon energy through the absorption coefficient (u,/p). The variation 
of this ratio with energy is shown in Fig. 5.1 (which is based on a value 
of 34.0 ev/ion pair for W,). It will be seen that between about 60 kev 
and 2 Mev the variations from a mean value of 3350 ergs/cm? per 
roentgen are not much greater than 10%. In this region of energies, 
therefore, the exposure rate provides a rough measure of intensity. 


roentgens (5.2) 


result from subjecting the material in question to an exposure of one roentgen.” 
Not too surprisingly, a certain amount of confusion resulted from this dual usage. 

Some clarification was achieved by referring in the second instance to a dose of one 
rep (roentgen-equivalent-physical), but a further difficulty arose when some workers 
based the rep on the energy deposited by one roentgen in air, and others based it 
on the energy deposited in soft tissue; the sizes of the unit in the two cases differ by 
about 10%. The introduction of the rad as a unit of absorbed dose removed this 
difficulty. 
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The theoretical relationship between exposure and quantity of radia- 
tion has been checked experimentally at a number of radiation energies 
(references 1-5 of Chapter 4). The method used in each case was to 
measure the absolute intensity at a point in the radiation beam with a 
calorimeter, and then to measure the exposure rate at the same point 
with a cavity chamber. (Cavity chambers will be discussed later in the 
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Fig. 5.1. Ratio of quantity of radiation to exposure dose as a function of photon 
energy. 


chapter.) In every case the measured number of ergs/cm? per roentgen 
agreed with the theoretical value within the experimental accuracy of 
2 to 3%. 

The free-air chamber. The definition of the roentgen requires that 
the ionization from all the secondary electrons originating in 1 em? of 
standard air be collected and measured. Obviously this cannot be done 
without at the same time collecting ionization from electrons originating 
in nearby volumes. However, one may make use of the equilibrium be- 
tween the x-rays and their secondary electrons to achieve the same end. 

Consider a parallel beam of x-rays of fixed cross-sectional area passing 
through air in the Z-direction with negligible attenuation (see Fig. 5.2). 
The intensity will be the same at all points along the beam. At a suffi- 
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cient distance from the source of radiation there will be an equilibrium 
in the Z-direction between the x-rays and the secondary electrons: as 
many electrons will stop in a layer Az as originate in that layer. Provided 
the layer extends laterally for a distance greater than the maximum 
range of the secondary electrons, the total ionization produced in it 
will be equal to the ionization produced by all the secondary electrons 
originating in the segment of the beam contained in Az; the ionization 
lost when electrons such as A (Fig. 5.2) leave the layer will just be com- 


YUsy B A 
ly 
fee 
cn ane 


M1111 1 1A 


Fig. 5.2. Principle of the free-air chamber. The total ionization produced in the 
layer Az by secondary electrons crossing it is equal to the total ionization produced 
by the electrons originating in Az. A, B, and C are typical secondary electrons. 
| 

pensated by the ionization gained from electrons such as B and C enter- 
ing the layer from outside. It follows that by collecting and measuring 
all the ionization produced in a layer of air of known thickness oriented 
at right angles to a beam of x-rays of known cross-sectional area, one can 
calculate the ionization resulting from the interactions in a unit volume 
of air and, hence, the exposure dose in roentgens. This principle is the 
basis of the standard free-air chamber, with which absolute measure- 
ments of exposure are carried out. 

The principal features of a typical free-air chamber are shown in 
- Fig. 5.3. The diaphragm, constructed of a heavy metal such as tungsten 
or gold, defines an x-ray beam of accurately known cross-sectional area, 
A. This beam passes between a pair of parallel plates in an air-filled 
enclosure. The upper plate is maintained at a high potential relative to 
the lower one, which is actually divided into three—a collecting electrode 
in the center and guard plates on either side. These guard plates help 
to keep the electric field perpendicular to the plates throughout the col- 
lecting volume, as do a set of guard wires or strips which are spaced at 
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equal intervals between the plates and maintained at a series of inter- 
mediate potentials. In the electric field arising from the potential dif- 
ference between the upper and lower plates all the ions of one sign pro- 
duced in the region between the dashed lines move to the collecting 
electrode. The value of the charge collected in this way is determined 
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Fig. 5.3. Principal features of a free-air chamber. 


to an accuracy of better than 0.1% by means of a sensitive measuring 
system. 

According to the argument presented above, the measured ionization, 
Q, is equal to that produced by the secondary electrons originating in 
the shaded volume (Fig. 5.3), provided the separation of the plates is 
great enough that no secondary electrons can strike them, and the dis- 
tance between the diaphragm and the collecting volume is greater than 
the maximum range of the secondary electrons. The shaded volume is 
proportional to the width L of the collecting electrode and to the area 
of the beam at the effective center of the volume. From the fact that 
the area of the beam varies as the square of the distance from the source 
of radiation, whereas the intensity varies inversely as the square of this 
distance, it can be shown that the charge produced by the secondary 
electrons generated per unit volume at the center of the diaphragm is 
very nearly equal to Q/AL. If the air density in the chamber is p g/cm®, 
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the exposure at the position of the diaphragm is therefore given by 


0.001293 
R= Q see 


AL ; roentgens (5.3) 
where Q is the charge in esu collected in the chamber, A is the area 
of the diaphragm in cm’, and L is the width of the collecting electrode 
in cm. 

In practice, the attenuation of the beam between the diaphragm and 
the collecting volume is not quite negligible, and a small correction must 
be made. Account must also be taken of the presence of water vapour 
in the air and the possible loss of ions by recombination. Finally, radia- 
tion scattered out of the main beam can cause a small amount of ex- 
traneous ionization in the collecting volume, and here again a small 
correction must be applied. When account is taken of all the relevant 
factors, a determination of exposure in roentgens with a free-air chamber 
should be accurate to better than 1% (2, 1). 

Most free-air chambers have plate separations of between 10 and 30 cm 
and are used in the conventional x-ray region of about 50 to 300 kv 
(3, 4, 5). At tube voltages below 50 kv considerably smaller chambers 
are generally used, and the electrode system is sometimes arranged in 
such a way that it can be operated at different distances from the dia- 
phragm (6), so that the correction for air attenuation can be measured 
directly. At very low tube voltages the extremely high intensities 
usually encountered make the complete collection of the ion current more 
dificult. Furthermore, the close correlation between exposure and 
absorbed dose in soft tissue which prevails above about 50 kev breaks 
down at these energies, so that the concept of exposure dose loses some of 
its usefulness. In this region of energies (sometimes called the grenz-ray 
region) it may prove more useful to measure intensities, either directly 
by means of a calorimeter or indirectly by means of some other form of 
totally absorbing detector (such as a scintillating crystal), and to derive 
the absorbed dose from the intensity. 

At high energies (above about 400 kv) difficulties arise because of the 
large ranges of the secondary electrons. Either one must make the 
chamber very large (7), in which case the problem of establishing a uni- 
form electric field strong enough for complete ion collection is very great, 
or one must operate the chamber at an air pressure of several atmos- 
pheres. The latter technique has been described in some detail by Taylor 
and Singer (8), who made measurements on the y-rays from radium with 
a chamber operated at 10 atmospheres; this technique also presents 
major difficulties, and has been used only rarely. 
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Free-air chambers are maintained in the standardizing laboratories of 
several countries. They constitute the national standards of exposure 
dose in the range from 10 to 300 kv. 

The cavity ionization chamber. As mentioned above, the use of 
free-air chambers at energies above about 400 kv becomes increasingly 
difficult because of the requirement that practically all the secondary 
electrons originating between the diaphragm and the collecting volume 
expend their full ranges within the air in the chamber. Evidently some 
modification in the method of measurement is necessary. 

One such modification leads to the principle of cavity ionization. 
Consider a large volume of air, uniformly irradiated by x-rays. At all 
points farther distant from the source of radiation than the maximum 
range of the secondary electrons, the secondary electrons will be in equi- 
librium with the x-rays, and the ionization produced in a small volume 
of air by all the electrons crossing it will be equal to the total ionization 
produced by the secondary electrons originating in that volume. The 
lonization per 0.001293 g of air under these conditions will be equal to 
the exposure in roentgens at the point of measurement: 


R = 0.001293J4 roentgens (5.4) 


where F# is the exposure in roentgens and J, 1s the ionization per unit 
mass of air, measured in esu/g. 

Now consider the effect of surrounding the small air volume from which 
the ionization is to be collected with a denser material having, apart 
from its density, the same physical properties as air (Fig. 5.4). The 





Fig. 5.4. Illustrating the effect of doubling the density of the medium surrounding 
a small air volume. The arrows represent secondary electrons. 


number of secondary electrons produced per cm? is proportional to the 
number of electrons present per cm? in the medium, i.e., to the electrons 
density. The specific energy loss of these secondary electrons is also pro- 
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portional to the electron density; their ranges are therefore inversely 
proportional to the electron density. Thus, if the density of the medium 
surrounding the air volume is doubled, secondary electrons will be able 
to reach this volume from only half as deep in the medium, but since 
twice as many will be produced per cm*, the number crossing the air 
volume will remain the same as before. It follows that the ionization 
in the air volume will be independent of the density of the material sur- 
rounding it,* and the large expanse of air necessary in the free-air cham- 
ber can be replaced by a very much smaller expanse of solid. 

A number of factors must be taken into account, however, before a 
measurement of the ionization in an air-filled cavity in a solid medium 
can be interpreted in terms of roentgens. These are discussed in the 
following sections. 

Equilibrium wall thickness. If the rate of ionization in a gas-filled 
cavity exposed to a beam of y-rays or high-energy x-rays is measured 
as a function of the thickness of solid material surrounding the cavity, 
the results can be represented by a curve of the sort shown in Fig. 5.5. 
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Fig. 5.5. Variation of the ionization in a gas-filled cavity with the thickness of the 
surrounding wall. 


The rate of ionization increases steeply as the thickness is increased from 
zero, then levels off to a broad maximum, and then decreases more 
gradually. The reason for the initial increase is that most of the ioniza- 
tion in the cavity is produced by secondary electrons originating in the 
walls; adding to the wall thickness results in more electrons crossing the 
cavity (layers 2 and 3 in Fig. 5.6), until a point is reached where electrons 
generated in the outermost layer (layer 4 in the figure) cannot reach the 
cavity. From this point on there is no further increase; the secondary 
electrons are in equilibrium with the primary radiation, and the ioniza- 


* A rigorous proof of this statement has been given by Fano (9). 
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wall of equilibrium thickness for 200 kv x-rays the attenuation is less than 
0.2%, whereas for Co y-rays it is a bit more than 1%. Some approxi- 
mate values of wall attenuation are given in Table 5.1. 

Stopping-power ratio. If one could surround the air volume in 
which the ionization is to be measured with a solid of exactly the same 
stopping power per electron as air, the observed ionization in the volume 
would be exactly the same as it would be in a similar volume in free air. 
Unfortunately there is no such “air-equivalent’”’ material, and a correc- 
tion for stopping power must always be made. Since the problem of 
stopping power affects measurements of absorbed dose by cavity ioniza- 
tion in much the same way as it affects measurements of exposure, it will 
be discussed here in some detail. 

The way in which the difference in stopping power between the air and 
the surrounding wall affects the observed ionization can be seen from the 
following rough considerations. Suppose an “air-equivalent’” wall 
around an air-filled cavity in a uniform x-ray field is replaced by a wall 
with a higher stopping power per electron. Since, as was pointed out 
above, the ionization in the cavity will not depend on the actual density 
of the material surrounding it, we may justifiably simplify the argument 
by choosing the same electron density for the new wall material as for 
the old. Then, if photoelectric interactions and pair production can be 
neglected, the radiation will produce the same number of secondary elec- 
trons per unit volume in the new wall. But because of the higher stop- 
ping power, these electrons will not travel as far in the wall material, and 
fewer of them will cross the air volume. The ionization will therefore be 
less by a factor Ry/Ra, where Ry is the range of an electron in the wall 
material, R4 is the corresponding range in the air-equivalent material, 
and the bar over the ratio indicates that a mean 1s taken over all ranges 
present. Since the distance an electron travels while its energy lies in 
any particular interval dE at E is inversely proportional to the stopping 
power at that energy,* the ratio Ry /R,, can be replaced by a correspond- 
ing ratio S4/Sy. Thus the ionization J4™ in the non-air-wall chamber is 
related to the ionization J44 in the air-wall chamber by an equation of 
the form 








Jan 





(The subscripts on the J’s refer to the gas, the superscripts to the wall 
material.) 

To obtain an exact expression for the mean stopping-power ratio, one 
must examine the situation in more detail. Consider a region of a solid 


* Strictly speaking, this is so only if the energy loss is continuous. 
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irradiated uniformly by x-rays. The energy per gram imparted to the 
medium at a point will be equal to the product of the flux of secondary 
electrons and the mass stopping power of the medium. Under condi- 
tions of “absolute” electron equilibrium this will be equal to the energy 
lost by the x-rays at that point, i.e., 


f $(E) mS dE = I,(ta/p) (5.6) 


where ¢(£) is the differential electron flux, ,,S is the mass stopping power 
(also a function of electron energy), J, is the spherical radiation intensity, 
and u,/p is the energy mass absorption coefficient for the radiation. The 
integral extends over all the possible energies of the secondary electrons. 

Now suppose a small air-filled cavity to be inserted at the point of 
interest. If the cavity is small enough so that most secondary electrons 
lose only small fractions of their energies in crossing it, the secondary- 
electron spectrum will not be appreciably altered by the presence of the 
cavity. In this case the ionization per unit mass of air in the cavity will 
be given by 


t 
Ja" = — f ou(E) mSa dE (5.7) 
A 


since the electron flux distribution will be that for the wall material, 
but the ionization under consideration is occurring in air. Wg, is the 
mean energy required to produce an lon pair in air. 

The ionization in an air-filled cavity in an air-equivalent wall will be 


t 
Pee 
Ja* = ; fou ma dE (5.8) 


The ratio of the ionization in an air-filled cavity in an arbitrary wall 
material to that in an air-filled cavity in an air-equivalent wall is, from 
Eqs. 5.7 and 5.8, 


7,4 = FE (5.9) 
A 

f b4(E) mSa dE ; 4(E) mSa dE 
and when the two cavities are exposed to the same radiation field, the 


integral in the denominator involving ¢4(E) can be replaced by one in- 
volving ¢y(E) with the aid of Eq. 5.6: 


(ua/ P)A 


(ua/ p)M 





f ta(E) mSa dE = f éu(E) mSy dE (5.10) 
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so that 


Ja# — (ta/p) f om (E) mS (mSa/mSu) dE 


Ja4 ~ (ua / p)A 





(5.11) 
f éu(E) mS OE 


The expression in square brackets is evidently an average of mS4/mSm, 
the ratio of the stopping powers in air and in the wall, weighted according 
to the energy loss over the electron distribution. In fact, comparison 
with Eq. 5.5 shows that it corresponds to the mean stopping-power ratio 
introduced earlier. Thus Eq. 5.11 can be written 





fa (ua/p)u BT Be) dat iad 
(ua/ P)A 
where 
ee f bu (E) mSu(mSa/mSar) dE 
nA/nSuo = —_____"__—_ (5.13) 
f om (E) mSu dE 


Absolute measurements of exposure by cavity ionization are of most 
importance at photon energies between 0.3 and 3 Mev. In this region 
the only significant x-ray interaction in light elements is the Compton 


effect, and 
Ha/p = Ga/p = eFa(Z/A)No (5.14) 


where .a, is the Compton absorption cross section per electron. In this 
case Eq. 5.12 reduces to 
au _ (4/A)u —— 


Ja (Z/A)a ————— (mS.a/mSu Fat 


= (Sa/oSu)Ja4 (5.15) 


where .S4/eSm is given by Eq. 5.13, with stopping powers per electron 
replacing mass stopping powers. Here the ionization depends only on 
the relative stopping powers per electron of the wall and air. 

From Eq. 5.4 (and the argument following it) and Eq. 5.12, the ex- 
posure dose in roentgens is related to the ionization in a cavity of ar- 
bitrary wall material under equilibrium conditions by 


R = 0.001293 44(1 + a) 
(ua/p)a Ja“ (1 + a) 


= 0.001293 ————_—— 
(ua/p)M (mS4/mSur) 





(5.16) 
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where the factor (1 + a) is inserted to represent the correction for at- 
tenuation in the wall. 

Evaluation of the stopping-power ratio. The calculation of the 
stopping-power ratio on the basis of Eq. 5.13 requires a knowledge of 
the relative (though not the absolute) distribution in energy of the 
secondary electrons crossing the cavity. Since the exact determination 
of this distribution is a fairly complex problem, a more approximate 
type of calculation, based on a number of simplifying assumptions, will 
be considered first. 

As can be seen from Eq. 5.13, the mean stopping-power ratio mS4/mSM 
is obtained by multiplying the value of the ratio at any particular energy, 
E, by the energy $(£) »Sx dE lost by electrons in an interval of energy 
dE at E, and averaging over all values of EZ. If one assumes that the 
electrons lose energy continuously in very small amounts, one can write 
for the average value of »Sa/mSm over the slowing down of a single 
electron of initial energy Eo 


ees 1 E 0 
(wSa/mSie)! = 5 J (mS.a/mSar) AE (5.17) 


where F£ is the energy of the electron at some point in the slowing-down 


process. Averaging in turn over the spectrum N(E,) of initial electron 


energies, one obtains 
E 


_ (mSa/mSau)’N (Eo) Eo dE 
mSA/mSa E 


N (Eo) Eo dEo 
0 


Em Ey 
[ "Ne f (mSt/mSu) dE dE 
fo 


Em 
N (Eo) Eo dEo 


(5.18) 


0 


Here E,, is the maximum energy with which secondary electrons are pro- 
duced. 

If the secondary electrons originate only in Compton interactions, 
N(Eo)dEp is just the relative number of Compton electrons with energies 
between Ey and Eg + dE produced by radiation of energy hy, as given 
by the Klein-Nishina formula. Values of N(Eo) for a wide range of 
photon energies are given in tabular form by Johns et al. (13), and in 
graphical form by Nelms (14). 
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The ratio mSa/mSm can be calculated as a function of energy with the 
aid of the stopping-power formula, Eq. 2.2. It should be noted that the 
ratio is influenced by the density effect, which will be significant in the 
solid material of the wall but not in the air in the cavity. As pointed out 
in Chapter 2, the effect can be taken into account in two parts: (1) the 
zero-energy effect by using values of the mean excitation potential J 
appropriate to the solid and the gas (Table 2.1), and (2) the finite-energy 
effect by subtracting the energy-dependent quantity 6 from the stopping 
power for the wall. 

Equation 5.18 has been used by a number of workers (15, 16, 17, 18) 
to calculate the mean stopping-power ratio for specific y-ray spectra. 
The author (17) has modified a simplified expression derived by Laurence 
(15) to take into account more recent information on stopping powers 
and the density effect. A number of stopping-power ratios for particular 
materials and y-ray energies, calculated with the aid of this expression, 
are presented in Table 5.2.* 


Table 5.2 
MEAN STOPPING-POWER RATIOS, mS4/mSmu 


hy (Mev) Graphite Polystyrene Lucite Bakelite Aluminum 


0.1 0.989 0.88 0.89 0.92 1.18 
0.2 0.991 0.89 0.89 0.93 1.16 
0.5 0.993 0.89 0.90 0.93 1.14 
1.0 0.996 0.90 0.90 0.93 1.13 
2.0 1.005 0.91 0.91 0.94 1.12 
3.0 1.014 0.92 0.92 0.95 1.12 


As will be shown in the next section, the numbers in Table 5.2 are 
based on a simplified view of cavity ionization, and can only be con- 
sidered approximate. Nevertheless, the more exact theory indicates 
that for cavities of the size commonly used in radiation work they should 
not be in error by more than a few tenths of a percent. Recent attempts 
to measure the stopping-power ratio under actual cavity-chamber condi- 
tions (19, 18) support this conclusion. 

Effects of 5-ray formation. According to Eq. 5.13, mSa/mSm can 
be interpreted as the ratio of the stopping powers in the air and in the 
wall material averaged over the electron spectrum crossing the cavity. 


* Attix has pointed out that Laurence’s paper contains an algebraic error which 
affects the values calculated in reference (17) toa small extent. However, comparison 
of certain of these results with values calculated more precisely by Attix and Ritz 
(18) shows the discrepancy to be less than half of one per cent. 
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The fact that a moving electron does not lose energy continuously in 
matter, but may transfer a considerable fraction of its energy to an 
atomic electron in a single encounter, must be taken into account in any 
detailed consideration of the stopping-power problem. 

The atomic electrons ejected with energies great enough so that they in 
turn can cause ionization are called 6-rays. In 1955 Spencer and Attix 
(20) and Burch (21) independently pointed out that the production of 
these 6-rays affects the mean stopping-power ratio in two ways. First, 
the spectrum of secondary electrons crossing the cavity 1s shifted toward 
low energies, since when a high-energy electron loses a large amount of 
energy to a 6-ray, the effect is to remove the high-energy electron from 
the spectrum and replace it with two electrons of lower energies. Sec- 
ondly, if an electron crossing the cavity loses energy to a 6-ray which 
then leaves the cavity, some fraction of the energy loss does not show up 
as lonization in the cavity. 

Spencer and Attix take these effects into account by defining an energy 
A, for which an electron can just traverse the gas in the cavity before 
stopping. This energy obviously depends on the size of the cavity and 
the density of the gas in it. Encounters involving transfers of energy 
greater than A are thought of as simply altering the electron spectrum; 
only losses of less than A are considered to be locally absorbed as ioniza- 
tion. All electrons with energies greater than A are considered to be part 
of the spectrum crossing the cavity, and the stopping powers are cal- 
culated to include only energy losses of A or less. The mean stopping- 
power ratio thus depends on A, and since A in turn depends on the size 
of the cavity, mS4/mSm also depends on cavity size, differing more from 
unity for smaller cavities. Spencer and Attix have calculated the varia- 
tion of the stopping-power ratio with cavity size for a number of radiation 
energies; their predictions are supported by the experimental results of 
Greening (22), Whyte (19), and Attix and Ritz (18). 

Burch’s theory (21) is still more elaborate. In addition to taking ac- 
count of the effects of 5-rays in a somewhat more detailed fashion, it 
makes provision for differences in the elastic scattering properties of the 
solid and the gas. (The greater the scattering in the gas, the longer will 
be the path of an electron crossing the cavity, and the greater will be the 
resulting ionization.) Unfortunately the evaluation of the scattering 
effect is extremely complicated, and it is difficult to obtain numerical 
results from the theory. 
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CHAPTER 


6 


Measurement 
of absorbed dose 


for x-rays 


It was pointed out in Chapter 1 that the quantity of significance in 
biological and medical work is not the amount of radiation passing 
through a point in matter, but rather the energy absorbed by the matter 
at that point. The measure of this energy absorption is the absorbed dose 
(defined in Chaper 1), and the unit is the “rad,’”’ equal to 100 erg/g. 
If Ey is the energy in ergs absorbed per gram of material, then the 
absorbed dose Dy is simply given by 


Du = 0.01Ey rads (6.1) 


The energy referred to is that deposited in matter at the point of in- 
terest, rather than that removed from the x-ray beam at that point. 
For example, when a high-energy x-ray undergoes a Compton interaction 
at a point in matter, an energy equal to that of the recoil electron is 
removed from the beam at that point, but the actual transfer of energy 
to the medium takes place all along the track of the electron. Thus x-ray 
interactions at a particular point contribute to the absorbed dose over a 
region near that point; similarly, the absorbed dose at a point is due to 
x-ray interactions in a region around the point. Only when the ranges 
of the secondary electrons are small compared to the volume of interest 
can one consider that the loss of energy from the radiation field and the 
absorption of this energy by matter occur at the same place. 

The absorbed dose is the quantity of primary importance in dosimetry. 
Such things as intensity and exposure are of secondary interest, or serve 
as a means of determining the absorbed dose. 

42 
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Absorbed dose by calorimetry. The most direct (and possibly the 
most difficult) method of measuring absorbed dose is by calorimetry. 
The basic idea is simple: to measure the amount of energy absorbed in a 
small mass of the material by detecting the attendant rise in temperature. 
Fig. 6.1 represents some of the basic features of such a measurement. 
A small volume of the medium at the place where the absorbed dose is to 
be measured must be isolated thermally from its surroundings, and the 
temperature rise in the volume must be measured with a sensitive device 


\ 
WN 


| i 
© > 


Fig. 6.1. Essentials of a calorimetric measurement of absorbed dose. 


such as a thermistor. A knowledge of the heat capacity of the medium 
allows the energy per gram corresponding to the rise in temperature, 
and hence the absorbed dose, to be calculated. 

In practice the experimental difficulties connected with such a meas- 
urement are great. Electrons escaping from the test volume must be 
exactly balanced by electrons entering from outside, which means that 
the material enclosing the test volume must not differ much from the ma- 
terial of the medium. The test volume must be small enough so that the 
radiation intensity does not vary appreciably over it. Unless the radia- 
tion intensity is high, the rate of change of the temperature is too small 
to be measured accurately. Finally, any part of the absorbed energy 
that goes to producing chemical changes in the medium does not appear 
as heat and hence is not detected. Only a very few measurements of this 
kind have been attempted (1, 2). 

Under certain conditions a different calorimetric technique can be 
used to obtain information about absorbed dose. If a uniform beam of 
radiation strikes a considerable area of the surface of an absorbing 
medium, the intensity near the central axis of the beam depends only on 
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the depth in the material, and not on the lateral position, i.e., if the z-axis 
is in the direction of the beam, the intensity is a function of z only. The 
decrease, d/,, in the plane intensity between the depth z and the depth 
z + dz is numerically equal to the energy lost by the beam in the layer 
dz. Then provided the ranges of the secondary electrons are small com- 
pared to the distance over which the intensity varies appreciably, the 
energy lost by the radiation in the layer dz can be equated to the energy 
actually absorbed by the matter in this layer, and the absorbed dose Dy 
can be obtained from the quantity d/,/dz: 


Duy = ——— rad 6.2 
iM oe (6.2) 


where p is the density of the material in g/cm’, J, is the plane intensity 
in erg/cm?-sec, and z is the depth in the medium in cm. 
On this basis the variation of absorbed dose with depth in any solid 
material can be measured with the arrangement shown in Fig. 6.2. The 
plane intensity is measured by means 
Radiation of a total-absorption calorimeter, as 
described in Chapter 4. By placing 
sheets of the material in question in 
front of the calorimeter, one can 





Absorbing ieee : ‘ 
4“ sheets measure the vanation of intensity 
PION with depth, from which d/,/dz can 
be obtained directly. 


A disadvantage of this type of 
measurement is that it does not take 
proper account of the effects of 
backscattered radiation. At energies 
where a significant fraction of the 
energy removed from the primary 
beam is backscattered the accuracy 
of the method is therefore low. The accuracy is also restricted by the 
fact that the result depends on the differentiation of an experimental 
curve. 

Absorbed dose by cavity ionization. A much more common 
method of measuring absorbed dose is by cavity ionization. Although 
in principle the method is less direct than the calorimetric one, in prac- 
tice the measurements involved are a good deal simpler. 

To begin with, the energy Ey absorbed per gram of material is related 
to the ionization per gram in the material, Jw“, through the mean 
energy, Wy, lost by electrons in the material per ion pair produced, 


Fig. 6.2. Arrangement for measur- 
ing variation of absorbed dose with 
depth in a solid medium. 
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(ergs) (esu) (ev) _y (ergs) 
a Ju! oe xX Wu ance xX 1.602 x 107 Gay: 
S 1 (ion pr.) 
4.803 X 107!°  (esu) 
= 3.336 X 10 VJa-! Wy ergs/g (6.3) 


From Eq. 6.1 it follows immediately that the absorbed dose in rads is 
related to the ionization per gram in the material by the equation 


Du = 3.336 X 107° ut Wa rads (6.4) 


Since it is not feasible to measure ionization in a solid or liquid, one 
must deal instead with the ionization in a gas-filled cavity in the medium. 
It follows from Eq. 5.7 that the ratio of the ionizations in the cavity gas 
and in the medium as given by 


up, J OE) mSa dE 
ee (6.5) 
W 
° [6sc(B) mS aE 


JG 
Ju 


where the subscript G is used to indicate that the gas in the cavity need 
not be air. Equation 6.5 can be written 


uw Wa | 6xe(B) mSue(wSa/mSue) GE 


Jg 
a 
z We f tu(E) mSu dE 
Ww ee 
= — (nSa/mSa) (6.6) 
We 
where 
f éu(E) mSit(mSa/mSu) dE 
ese = (6.7) 
f éu(E) mS dE 


The mean stopping-power ratio »Sg/mSm is identical in form to that 
introduced in connection with the measurement of exposure by cavity 
ionization and defined by Eq. 5.13. Equation 6.6 simply expresses the 
fact that since the same electron flux traverses the gas in the cavity and 
the medium around it, the resulting ionizations per unit mass are directly 
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proportional to the mass stopping powers and inversely proportional 
to the mean energies expended per ion pair produced. Combining Eas. 
6.6 and 6.4 gives the absorbed dose in the medium 


Du = 3.336 X 107° (mSca/mSm) JG! We rads (6.8) 


The quantity (mS¢/mSmu)' is the same as the quantity s,, defined in 
Section 5 of the 1956 Report of the ICRU (3). Tables 3, 4 and 5 of that 
report present data relevant to the computation of mass stopping-power 
ratios when the gas involved is air. In general, the quantity can be cal- 
culated by the methods outlined in Chapter 5 and illustrated by the re- 
sults in Table 5.2. The application of Eq. 6.8 also requires a knowledge 
of Wg. Values of this quantity for a number of gases are listed in Table 
2.4 (Chapter 2). 

Equation 6.8 is one expression of the well-known Bragg-Gray cavity 
principle. This principle, which was first enunciated by Bragg (4) and 
later proposed independently by Gray (5, 6), states that the energy given 
up by secondary electrons to unit volume of a solid medium is equal to 
the product of (a) the ionization per unit volume in a small gas-filled 
cavity in the medium, (b) the mean energy expended by the electrons in 
producing an ion pair in the gas, and (c) the ratio of the specific energy 
losses of the secondary electrons in the medium and the gas. This result 
follows directly from the assumption that the secondary electron dis- 
tribution at the position of the cavity is unaffected by the presence of 
the cavity (and that the energy losses are continuous). 

It should be noted that in the measurement of absorbed dose the re- 
quirement of equilibrium between the radiation and the secondary elec- 
trons, so important in the measurement of exposure dose, no longer holds, 
since both the absorbed dose and the ionization from which it is cal- 
culated are determined by the flux of secondary electrons at the point of 
interest. At a point near the surface of the medium where the secondary 
electrons have not yet reached equilibrium with the incoming radiation, 
the ionization will be less than at a point deeper in the medium, but so, 
of course, will the absorbed dose. The ionization measurement will still 
be valid. 

A less obvious matter may also bear mentioning. It was stated earlier 
that one condition for a valid measurement is that the cavity be small 
enough so that the electron distribution crossing it will not be appreciably 
different from that in the adjoining medium. This condition becomes 
less restrictive when the walls and gas are closely similar in composition. 
In this case a measurement of cavity ionization, even in a cavity of con- 
siderable size, correctly yields the average absorbed dose in the material 
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immediately adjacent to the cavity. However, because of the attenua- 
tion of the radiation in the medium, this may not be quite the same as 
the absorbed dose that would have been delivered to the medium at the 
position of the center of the cavity if the cavity had not been there: if one 
imagined the ionization to be measured in progressively smaller cavities 
centered about the same point, one would expect that the average flux 
of electrons crossing the cavity would be progressively smaller because of 
increased attenuation, and the ionization per gram of gas in the cavity 
would be progressively less. Thus if the cavity in which the actual 
measurement is made is large enough, it may be necessary to reduce the 
observed ionization somewhat to correct for this effect. 

The same considerations apply to the measurement of the exposure at 
@ point in matter by cavity ionization: the exposure at the center of the 
cavity will tend to be slightly greater than it would be at the same point 
if the cavity were not there. 

Tissue-equivalent case. Attempts have been made to measure the 
absorbed dose in tissue through the use of a cavity chamber having a 
tissue-equivalent wall and filled with a tissue-equivalent gas. The 
absorbed dose in the tissue Dr is thus equal to the absorbed dose in the 
wall material Dy. The stopping powers in the gas and the solid differ 
only by virtue of the density effect in the solid. At radiation energies 
below about 1.5 Mev this difference is not likely to be greater than 2%; 
to an accuracy of this order one can therefore write 


Dr = Du = 3.336 X 105(.wSG/mSu) WG“ We 
= 3.336 X 10 °Jg” Weg rads (6.9) 


and the only parameter that must be known independently is the quan- 
tity Wg for the tissue-equivalent gas. 

Cavity with arbitrary wall. The tissue-equivalent case is the 
simplest theoretically, although by no means the simplest to realize 
experimentally. At the other extreme is the case where an air-filled 
cavity of some arbitrary wall material is imbedded in a medium such as 
tissue in order to determine the absorbed dose in the tissue. In this case 
the absorbed dose Dy in the wall material is given by Eq. 6.8. Since the 
radiation energy absorbed per gram of material is proportional to the 
energy mass absorption coefficient ua/p, the absorbed dose Dr in the 
surrounding tissue is given by 





(ua/ p)T 


(Wa/p)ae 





Dr = Dy (6.10) 


18 PRINCIPLES OF RADIATION DOSIMETRY __—“C—sSCS—S 


The bars over the absorption coefficient indicate averages over the 
radiation spectrum, i.e., 





J (oa/eyTa(ho) a) 
fp eee eae (6.11) 
f I, (hy) d(hv) 


where /,(hv) is the spectral distribution of the spherical intensity. 
In the energy region where only Compton interactions are significant, 
the energy mass absorption coefficient is given by Eq. 5.14 and 


(ue/p)r  (4/A)r 











= (6.12) 
(ua/p)m = (Z/A)u 
for any photon energy. In this case 
(Z/A)r 
= 6.13 
T MT Ayu (6.13) 


This last equation simply states that when the absorbed dose arises 
solely from Compton interactions, which occur in direct proportion to 
the number of electrons present, the absorbed dose in a particular ma- 
terial is proportional to the number of electrons per gram in that ma- 
terial, and hence to Z/A. 

For a measurement of the sort described above to be valid, the measur- 
ing chamber must be small, and it must not disturb the radiation in- 
tensity at the point of measurement. These conditions are most readily 
met with a condenser chamber, consisting usually of a hollow cylinder 
of a low-atomic-number conducting material enclosing a small volume 
of air, and an inner electrode insulated from the outer shell. The cham- 
ber behaves like a small condenser, and before the measurement it is 
charged up by the application of an appropriate voltage. When it is 
placed in the medium and exposed to the radiation, the ionization in the 
air volume conducts charge from one electrode to the other. The drop 
in voltage occurring during the period of irradiation multiplied by the 
capacitance of the chamber, gives the charge released in the air by the 
ionization. The absorbed dose in the wall can then be obtained from 
Eq. 6.8, and in the surrounding tissue (or other material) from Eq. 6.10 
or Eq. 6.13. 

It should be noted that in order for Eq. 6.8 to give the absorbed dose 
in the wall, all the electrons crossing the air volume must originate in 
the wall material and none in the surrounding medium. This means 
that the wall must be thick enough to prevent electrons from outside 
the chamber reaching the air volume. If the wall is not thick enough, 


MEASUREMENT OF ABSORBED DOSE FOR X-RAYS 79 





the observed ionization lies between that for a cavity in the wall material 
and that for a cavity in the surrounding medium; for an extremely thin 
wall the ionization corresponds wholly to the medium. 

Comparison between ionization and calorimetric methods. 
A useful check on the absolute methods of determining absorbed dose 
discussed above is provided by a comparison of calorimetric and ioniza- 
tion measurements made at the same point in the radiation field. In one 
such comparison (2) the absorbed dose in graphite at a point in a beam of 
Co® +-rays was determined first from the rise in temperature of a hollow 
graphite cylinder at that point, and then from the ionization in an air- 
filled cavity in a similar graphite cylinder at the same point. The two 
measurements were found to agree if a value of 33.0 ev/ion pair was 
chosen for W4. This agrees within 3% with the value of 34.0 obtained 
by a number of other methods. Whether the difference is significant 
will have to be decided by further experiment. 

The extrapolation chamber. Under certain conditions it is useful 
to be able to deduce the ionization per gram of gas that would occur in 
a cavity of vanishingly small volume. This is so, for example, when the 
ranges of the secondary electrons are very short, or when measurements 
are required in a region where the electron flux changes rapidly with 
depth in the material, as it does near the surface of an irradiated medium. 

One solution to this problem was devised by Failla (7), who developed 
a flat, disc-shaped chamber of adjustable electrode separation, called an 
extrapolation chamber. The essential features of such a chamber are 
indicated in Fig. 6.3. The collecting electrode is a thin disk of conducting 
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Fig. 6.3. Extrapolation chamber. 


material on a thick base of a material that simulates the medium of in- 
terest. Separated from the collecting electrode by a non-conducting 
gap is a guard electrode, which is maintained at or near the potential of 
the collector. The high-voltage electrode, generally in the form of a thin 
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foil, can be moved up or down with respect to the collector by means of a 
screw arrangement. The ionization current is measured as a function of 
the distance between the electrodes, and the ionization per unit mass of 
gas corresponding to a vanishingly small gap can then be obtained by 
extrapolation. (The validity of this procedure when the measurements 
are carried to very small gaps has recently been brought into question by 
the work of Spencer and Attix and of Burch referred to in the previous 
chapter. ) 

In addition to their use in x- and y-ray dosimetry (8), extrapolation 
chambers have been applied to the measurement of absorbed dose rates 
from §-ray sources (9, 10), and from sources uniformly distributed 
through matter (7). 

Calculation of absorbed dose from exposure dose. In principle 
the direct determination of absorbed dose always involves an absolute 
measurement at the point of interest in the medium. It is not possible 
to have a secondary detector calibrated to read directly in rads, since the 
absorbed dose depends on the properties of the particular medium in 
which it is measured.* On the other hand, cavity chambers are readily 
calibrated to read exposure dose in roentgens, and can be made to do so 
independently of the radiation spectrum over a considerable range of 
energies. Furthermore, extensive tables have been compiled giving the 
exposure as a function of depth in a tissue-like material for a given ex- 
posure in air at the position of the surface. It is therefore often con- 
venient to determine the exposure experimentally and then calculate 
the corresponding absorbed dose at the place of interest. Indeed, the 
fact that this conversion from exposure dose to absorbed dose can be 
carried out with reasonable accuracy over a wide range of energies—and 
that the conversion factor for soft tissue is very nearly independent of 
energy over much of this range—is one of the chief reasons for the use 
of the concept of exposure dose. 

The problem then is: given the exposure at a point in matter of known 
composition, to calculate the absorbed dose at that point. It is necessary 
for the calculation that there be electron equilibrium and that the radia- 
tion intensity not vary appreciably over distances of the order of the 
ranges of the secondary electrons. 

In Eq. 5.2 the exposure #& was related to the radiation intensity J,. 
Under the conditions specified above, the absorbed dose Dy, can also be 
related to the intensity: 


Dy = 0.01 FE x = 0.01 (ua/p)arl ot rads (6.14) 


* It is possible, of course, to calibrate a detector to read rads in one particular ma- 
terial, e.g., ‘‘rads in tissue” or “‘rads in bone.”’ 
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Combining Eas. 5.2 and 6.14 gives 
(Ha / p)M 
Ba/ p)A 


and if W, is put equal to the currently accepted value of 34.0 ev/ion 
pair, this reduces to 


Dy = 0.02580 





Wak rads (6.15) 


(ua/ p)M 


Du = 0.877 
(ua/p)A 





R rads (6.16) 


Thus, under the conditions mentioned above, the conversion factor from 
exposure in roentgens to absorbed dose in rads depends only on the ratio 
of the energy absorption coefficients of the medium and air. Where 
radiation of more than one photon energy is present, the absorption co- 
efficients must, of course, be averaged over the radiation spectrum. 

Between about 0.1 Mev and 3 Mev the only significant x-ray inter- 
action in light elements is the Compton effect, and the absorption co- 
efficient ue/p can be replaced by the Compton absorption coefficient 
o,/p. The latter is directly proportional to the number of electrons per 
gram in the medium, and hence to the mean value of Z/A. In this region 
of energies Eq. 6.15 can therefore be replaced by 


(2/A)u 
(Z/A)a 


where Z/A = 2f,(Z/A):, f; being the fraction by weight of the ith ele- 
ment in the medium. Some values of Z/A and the corresponding values 
of the conversion factor from roentgens to rads (based on Wy, = 34.0 
ev/ion pair) are listed in Table 6.1. The compositions assumed for 
the various materials are those listed on page 12 of reference 3. 





Du = 0.02580 WaR rads (6.17) 





Table 6.1 
CONVERSION Factors FROM ROENTGENS TO Raps IN SEVERAL MATERIALS 
0.1 to 3.0 Mev 
. as (Z/A)m 
1 : ——= 

Materia Z/A 0.87; (Z/A)a (rads/r) 
Air 0.499 0.877 
Water 0.555 0.97, 
Muscle 0.550 0.96, 
Bone 0.530 0.93; 


As stated above, the values of the conversion factors are practically 
constant between 0.1 and 3 Mev. At photon energies below about 0.1 
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Mev the photoelectric effect results in a variation of the conversion factor 
with energy; at high energies pair production leads to a similar variation. 
These variations are particularly pronounced in bone, which contains 
over 10% by weight of calcium (Z = 20). The variation of the absorbed 
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Fig. 6.4. Absorbed dose in rads per roentgen of exposure as a function of photon 
energy. 


dose per roentgen with energy for a number of materials of medical in- 
terest is shown in Fig. 6.4. The values for a number of typical radiation 
spectra, taken from reference 3, are listed in Table 6.2. 

The conversion factors derived above enable one to calculate the 
absorbed dose at a point in matter from the exposure at that point. 
Under most standard conditions, however, it is not necessary actually 
to measure the exposure inside the medium. Extensive measurements 
have been carried out, relating the exposure at various depths in a tissue- 
like medium to the exposure measured at the position of the surface in 
the absence of the medium (sometimes called the ‘‘air dose’). Such 
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Table 6.2 


CONVERSION Facrors FROM ROENTGENS TO Raps 
FOR TYPICAL RapIATION SPECTRA 


rads/r 
Filter HVL 

Radiation (mm) (mm) Water Muscle’ Bone 
100 kv x-rays 0.18 Cu 5.5 Al 0.91 0.94 3.1 
150 0.075 Cu 0.2 Cu 0.92 0.94 2.7 
200 0.20 Cu 0.5 Cu 0.94 0.95 2.0 
250 0.17 Cu+ 3.0 Al 1.0 Cu 0.95 0.95 1.76 
250 0.9 Cu + 3.0 Al 2.0 Cu 0.96 0.96 1.42 
280 — 2.5 Cu 0.97 0.96 1.22 
400 — 4.2 Cu 0.97 0.965 1.11 
Cs!87 y-rays — — 0.97, 0.96, 0.93 

(0.66 Mev) 
Co y-rays — — 0.975 0.96; 0.93 

(1.25 Mev) 


depth-dose * data have been compiled for a large number of radiation 
spectra, beam areas, and distances from the source of radiation to the 
surface of the medium (11, 12). A few representative depth-dose curves 
for x-rays are shown in Fig. 6.5. 

Under radiation conditions for which depth-dose data are available, 
the simplest procedure for determining absorbed dose at a point in tissue 
is to measure the exposure rate in air at the point where the surface of 
the tissue will be, determine from the depth-dose data the resulting ex- 
posure rate at the required depth in tissue, and then convert this exposure 
rate to the corresponding absorbed-dose rate in tissue with the aid of 
Eq. 6.16 or Table 6.2. The absorbed dose rate in fat or bone at this point 
can also be obtained simply by using the conversion factors appropriate 


* The use of the term ‘‘dose”’ in this connection dates back to the time when the 
distinction between the concepts of exposure and absorbed dose was not always ex- 
plicitly made. The standard measurements have all been made with air-filled cavity 
chambers, and, where the results are expressed as a percentage of the number of 
roentgens measured in air at the surface, they must be interpreted as giving the ex- 
posure in roentgens. However, at the higher energies, where a thin-walled chamber 
has been used and the build-up of the ionization near the surface has been detected, 
the readings must be interpreted as giving relative values of the absorbed dose in 
tissue. To use these in conjunction with a measurement of exposure in air, one must 
still deduce from the latter the exposure at, say, the maximum of the depth dose curve 
and then calculate the corresponding absorbed dose, thereby “normalizing’’ the curve 
at that point. 
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Fig. 6.5. Depth-dose curves for x-radiation in water. (Distance from source of 
radiation to surface of water, 80 cm; area of beam at surface of water, 100 cm?.) 


to these substances. It should be noted that in a small volume of bone 
or fat embedded in tissue, the exposure will differ very little from that in 
straight tissue at the same depth, but the absorbed dose may differ con- 
siderably from that in tissue. 
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CHAPTER 


d 


Dosimetry 
of charged particles 


The techniques of charged-particle dosimetry differ from those of 
x-ray dosimetry partly for historical reasons, but mainly because of 
differences in the rates of energy loss. Except at very low energies, an 
x-ray beam loses only a few per cent of its energy in traversing a centi- 
meter of light material: the absorbed dose can be considered essentially 
constant over the volume of a detector of reasonable size. On the other 
hand a typical beam of charged particles loses a large fraction of its 
energy in a centimeter of matter, and hence the absorbed dose varies 
greatly in a short distance. This fact has led to the development of 
special methods for charged-particle dosimetry. 


ELECTRONS 


Sources. High-energy electrons are obtained in two principal ways: 
as electron beams from particle accelerators, and as 6-rays from radio- 
active nuclides. Particle accelerators include Cockcroft-Walton ma- 
chines, Van de Graaff generators, linear accelerators, betatrons and 
synchrotrons. The first two of these produce electrons with energies 
up to a few Mev, linear accelerators commonly operate between about 
4 and 20 Mev, and betatrons and synchrotrons between 10 and 300 Mev. 
All of these machines produce well defined beams of nearly monoenergetic 
electrons. Beta-ray sources, on the other hand, emit electrons in all di- 
rections and with a continuous distribution of energies. Beta emitters 
with end-point energies above 3 Mev are rare, and the great majority 
lie below 2 Mev. Since the range of 2-Mev electrons in water is only 
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about 1 cm, it is evident that the effects of discrete 6-ray sources are con- 
fined to small thicknesses of matter. 

Intensity measurements. The intensity of a beam of electrons is 
determined by the flux (electrons/cm?-sec) in the beam and the energy 
of the electrons. The electron energy is frequently known from the con- 
ditions of operation of the accelerator, or it can be determined by 
magnetic-deflection or range measurements. The beam flux is deter- 
mined either with a Faraday cup or by ionization in air. 

The essential features of a Faraday cup are shown in Fig. 7.1. The 
electrons enter through a thin window and strike a metal block which 
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Fig. 7.1. Schematic diagram of a Faraday cage. 


is supported on insulators in an evacuated chamber. The block is thick 
enough to stop the incident electrons completely, and hence it acquires 
one negative charge for every electron passing through the window. 
(It may be necessary to correct for electrons scattered back out of the 
block before being stopped.) The current flowing from the block is thus 
equal to the current in the portion of the beam passing through the 
window. When combined with a knowledge of the beam (or window) 
area, and electron energy, this measurement yields the beam intensity. 

In theory an ionization chamber can also be used for an absolute 
determination of the intensity in an electron beam. If a beam of ¢ 
electrons/cm?-sec of energy E passes through a layer of gas in which the 
specific energy loss of the individual electrons is —dE/dz and the mean 
energy to produce an ion pair is W, the ionization produced per unit 
volume in a time ¢ should be given by 


_ g¢t(—dE/dz) 


J 
7 W 


(7.1) 
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If the electron energy is known, the specific energy loss can be computed 
from the stopping-power formula; a measurement of the ionization 
should thus yield a value for the electron flux, ¢, and the intensity is 
simply I = ¢E. [The quantity of radiation (erg/cm?) in an electron 
beam corresponding to an ionization density of 1 esu/cm® in air has been 
calculated as a function of electron energy by Laughlin (1).] 

The ionization density Jy can be measured in a parallel-plate ioniza- 
tion chamber similar in construction to the free-air chambers used with 
x-rays. In practice scattering of the electrons makes necessary rather 
large corrections, and when high accuracy is required it may be necessary 
to calibrate the chamber against a Faraday cup. In this case the ioniza- 
tion chamber becomes a secondary instrument. 

A cavity chamber with walls thick enough to establish equilibrium 
between the primary electrons and their 6-rays can also provide a means 
of measuring the ionization density. Again the accuracy is improved by 
calibration against a Faraday cup. 

The intensity measurements described above are applicable to directed 
beams of high-energy electrons. The intensity due to a point source of 
B-rays is much more difficult to measure because it is such a strongly 
varying function of position. Beta-ray sources distributed over a plane 
area are sometimes used (in eye applicators, for example) and in this case 
the intensity can be determined with the aid of an ionization chamber 
with a shallow, disk-shaped volume parallel to the plane of the source (2). 
Since the energy spectrum and angular distribution of the 8-rays crossing 
the chamber volume are not accurately known, the intensity cannot be 
inferred from a single ionization measurement at the surface. Instead 
the ionization occurring in a thin layer of air is measured at various 
depths in a suitable absorber (which constitutes the upper and lower 
plates of the chamber), the corresponding rate of energy dissipation per 
unit thickness of absorber is calculated for each depth, and the results 
are integrated from zero to a depth equal to the maximum range of the 
B-particles. This integral gives the total energy per second deposited by 
the 8-particles in the absorbing medium, which must be equal to the 
total energy per second emitted by the source. When this is divided by 
the source area, the result is the average plane intensity at the surface 
of the source. 

Absorbed-dose measurements. The most direct way to measure 
the absorbed dose due to a flux of electrons would be to determine calo- 
rimetrically the heat developed per unit mass of material at the point 
of interest, but as yet this has not proved experimentally feasible. A 
second method is to measure the particle flux as described above and 
then use the stopping-power formula to calculate the energy absorbed 
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per gram of the material of interest. This procedure is best suited to 
thin layers of material, for which the beam energy is well defined; in 
thick layers, where the primary electrons are scattered and degraded in 
energy, the calculation is less straightforward (3). As in the case of 
x-rays, the absorbed dose is most commonly measured by cavity ioniza- 
tion. 

The relationship between the absorbed dose at a point in a medium 
traversed by a flux of electrons and the ionization in a small gas-filled 
cavity at that point is the same as for a medium exposed to x-rays, since 
in both cases the absorbed dose is determined solely by the electron flux 
and energy distribution. Thus the dose in a medium M is given in terms 
of the ionization in a cavity filled with a gas G by Eq. 6.8. Again the 
quantity mSc/mSm is the ratio of the mass stopping power of the gas to 
that for the medium, averaged over the electron spectrum crossing the 
cavity. One of the chief sources of error in such a determination lies in 
the evaluation of this stopping-power ratio. This is particularly so at 
primary energies above a few Mev, where the presence of the density 
effect in the solid makes the ratio increasingly energy-dependent (4), and 
hence increasingly sensitive to the exact energy distribution assumed 
for the electrons crossing the cavity. Here the calculations of Spencer 
and Fano (5) on energy distributions of electrons slowing down in matter 
have been helpful. 

The way in which the absorbed dose due to an electron beam varies 
with depth in matter is often of interest. Such depth-dose curves are 
usually determined with a small cavity chamber, the ionization current 
being assumed proportional to the absorbed dose. Typical depth dose 
curves for a few beam energies (6, 7, 8) are shown in Fig. 7.2. The 
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Fig. 7.2. Depth-dose curves for electrons in water. 
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initial increase in the curves is due partly to the build-up of secondary 
(knock-on) electrons (which are as effective in producing ionization as 
the primaries), and partly to scattering of primary electrons (see p. 4). 
The exact shape of the depth dose curve at a particular energy depends 
somewhat on the area of the beam. The factors governing the penetra- 
tion of electrons in matter have been elucidated in theoretical studies 
by Spencer (3) and Adawi (9). 

The dosimetry of 8-rays has followed somewhat different lines. Beta- 
ray sources tend to be distributed over an area or throughout a volume. 
A source of any size and shape can be considered as made up of an ap- 
propriate array of elementary point sources. If the dose distribution 
around a point source is known, the dose at any position due to an actual 
source can be deduced by summing the contributions of the component 
point sources. The problem is thus reduced to finding the dose distribu- 
tion around a point source of $-rays in matter. 

Since the integrations involved in calculating the dose distribution 
from an extended source are greatly simplified if the point-source func- 
tion has a simple analytical form, a number of approximate but mathe- 
matically convenient expressions have been used for such calculations. 
For example, since the attenuation of 8-rays is roughly exponential over 
a considerable portion of the range, and since, apart from attenuation, 
the flux falls off as the square of the distance, a function of the form 


D(r) = ke~#*/r? (7.2) 


has been used (10, 11). Here r is the distance from the source and u is a 
rough empirical absorption coefficient dependent on the energy of the 
B-particles and on the absorbing material. 

The determination of a more exact point-source function is made 
difficult by the extremely non-uniform nature of the dose distribution 
around a point source of 6-rays. In a solid or liquid the dose falls off to 
zero within a few millimeters, and consequently the dimensions of the 
measuring instrument must be prohibitively small. One way to get 
around this difficulty is to use an extended plane source and make meas- 
urements with a parallel-plate extrapolation chamber (13). Since in this 
case the dose is constant over any plane parallel to the source, it is 
possible to use a chamber of reasonably large area, while the extrapola- 
tion technique allows results to be obtained for vanishingly small cavity 
thicknesses. It is thus possible to obtain a fairly accurate plane-source 
dose function, and it is then only necessary to deduce what form of point- 
source function this corresponds to. The result obtained by Loevinger 
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(12) for high-energy 8-rays such as those of P®* was 
D(r) = k/(ur)? r<l/u 

D(r) = ke!—#"/ur r>1/p (7.3) 


This expression falls off much less steeply at small distances than the 
simple exponential expression of Eq. 7.2, in agreement with experiment. 
The value of the constant uw for a particular source and material is not 
necessarily the same as in Eq. 7.2. 

Another method of determining the point-source function is to measure 
the dose distribution around a source in air. In this case the absorption 
takes place over distances of the order of meters rather than millimeters, 
and ionization chambers of reasonable size can be used. An analysis of a 
variety of results obtained by this method (13) shows that they can be 
fitted by a function that agrees with Eq. 7.3 for high-energy sources and 
takes a slightly more complicated form at lower energies. Using this 
function Loevinger, Japha and Brownell (14) have made a detailed 
analysis of dose distributions for a variety of source configurations. 

It has recently been demonstrated (15) that the attenuation curves 
obtained by interposing plane absorbers between a 8-ray source and a 
counter can be accurately represented by a binomial function of the form. 


N/No = (1 — 2/R)" (7.4) 


where N is the counting rate corresponding to a thickness x of absorber, 
No is the rate for no absorber, RF is the maximum range of the 8-particles, 
and n is a constant of the order of 3. This function has the advantage 
that it falls to zero at a thickness equal to the range, in contrast to the 
functions involving exponential factors which extend to infinite thick- 
nesses. However, the simple binomial expression of Eq. 7.4 fails to give 
an accurate representation of the point-source dose function, particularly 
at small distances. 


HEAVY CHARGED PARTICLES 


Heavy charged particles have found only limited use in situations re- 
quiring dose measurements. This is: because at the energies readily 
available the specific energy loss of a heavy particle is so high that it can 
penetrate only a very small thickness of matter—about 1 mm of water 
for a 10-Mev proton and 8 cm for a 100-Mev proton. The full energy 
of a beam of particles thus appears in a small volume, and hence heavy 
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particles are only useful where a large absorbed dose is required in a very 
limited region. 

There is another difference between the behavior of heavy particles 
on the one hand and x-rays and electrons on the other. The biological 
effect of a given absorbed dose of radiation depends on the average LET 
of the charged particles transferring the energy to the medium. The 
effectiveness of a particular kind of radiation in producing a particular 
biological effect is characterized by a number called the RBE (relative 
biological effectiveness), which is inversely proportional to the number 
of rads required to produce the effect. An RBE of unity is generally as- 
signed to 200 kvp x-rays; on this basis the RBE of another type of radia- 
tion (for a particular biological effect) is equal to the number of rads 
that must be delivered by 200 kvp x-rays to produce the effect, divided 
by the number of rads delivered by the other radiation to produce the 
same effect. For most types of biological effect the RBE of radium and 
other high-energy y-rays is somewhat less than unity, 0.7 being a typical 
value. For heavy charged particles, which have much higher LET 
values, the RBE may be very much greater than unity, and in fact 
values between 3 and 20 arecommon. This is also true of fast neutrons, 
which transfer energy to matter mainly through recoil nuclei. 

For certain purposes it is convenient to define an RBE dose, equal to 
the absorbed dose in rads multiplied by the RBE. The unit of RBE dose 
is the rem. Thus 


RBE dose (rems) = absorbed dose (rads) X RBE 


Since the biological effect of a given dose may depend on other factors 
besides the LET, RBE and RBE dose can only be considered as some- 
what crude concepts. They have been found useful for such applications 
as specifying permissible doses for radiation workers. 

Sources. Radioactive substances do not emit protons, and the 
a-particles emitted in radioactive transformations have such short ranges 
that they are of little biological interest, except where the a-emitter is 
actually incorporated into living material. The principal sources of 
heavy charged particles are therefore high-energy accelerators. Among 
these are the Cockcroft-Walton and Van de Graaff machines for energies 
up to a few Mev, and cyclotrons and linear accelerators for higher ener- 
gies. All these machines produce directed beams of heavy particles— 
most commonly protons, although deuterons and a-particles can also 
be accelerated. In the first two machines the beam current is usually 
continuous; in the last two it 1s pulsed. 

Intensity and absorbed dose. The methods of measuring intensity 
and absorbed dose for heavy-particle beams are largely the same as for 
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electron beams. Such differences as there are arise mainly from the fact 
that scattering and bremsstrahlung are much less important for heavy 
particles, whereas absorption in any material between the source and the 
detector is much more important. Intensity measurements are made 
with a Faraday cup similar to the kind used for electron beams, or with 
a, parallel-plate ionization chamber. Ionization-chamber measurements 
can be made to yield values of the absorbed dose: the beam is allowed to 
pass through a suitable gas, the ionization from a known volume is meas- 
ured, and from this the energy loss that would take place per gram of 
tissue is calculated. Ionization chambers for measurements of this kind 
must be designed so that the beam passes through a minimum of absorb- 
ing material in traversing the chamber. 

At energies above about 20 Mev another method of measuring particle 
flux—and, hence, indirectly, absorbed dose—becomes practicable for 
heavy particles. This method takes advantage of the fact that at these 
energies heavy particles can initiate nuclear reactions that result in the 
production of radioactive nuclides. If the number of radioactive nuclei 
resulting from an exposure to the beam can be measured by absolute 
counting methods and if the activation cross section is known, the num- 
ber of particles arriving per cm? per second in the beam can be de- 
termined from the equation 


S = $ee(t — 7.5 
yo eee (7.5) 


where A is the number of radioactive atoms at the end of an irradiation 
of duration ¢, N is the total number of atoms exposed to the beam, ¢ is 
the average flux in particles per cm? per sec in the beam, and ) is the 
decay constant of the radioactive nuclide. The method does not yield a 
high absolute accuracy, but it has been used successfully to check the 
reliability of the ionization method (16). 

The depth-dose curve for a charged particle beam is similar to the 
Bragg curve for particles of that energy (see Chapter 2). It is dis- 
tinguished by an increase toward the end of the range, corresponding 
to the increase in the specific energy loss with decreasing energy. 
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Neutron 
dosimetry 


The deposition of energy in matter by fast neutrons bears some 
resemblance to that by x-rays in that the loss of energy occurs in 
relatively large but infrequent transfers to secondary charged particles, 
which then deposit the energy in the medium by ionization and excita- 
tion. Neutrons, like x-rays, are thus able to carry energy to consider- 
able depths before transferring it to matter. On the other hand, the 
secondary particles that actually deposit the energy in the medium are 
heavy charged particles rather than electrons, and hence fast neutrons 
resemble heavy charged particles rather than x-rays in their biological 
effects, being characterized by RBE values of the order of 10. 

Sources. The neutron is not a stable particle in the free state. If it 
avoids capture by a nucleus, it will decay spontaneously into a proton 
and an electron, the half-life of the process being 13 minutes. Any 
neutron source must therefore be based on the production of neutrons 
in a nuclear reaction. Furthermore, since neutrons are electrically 
neutral, they cannot be accelerated to higher energies like charged 
particles; their energies are governed solely by the energy balance of 
the reaction in which they are produced. 

One of the most commonly used neutron-producing reactions is that 
resulting from the interaction of a-particles with beryllium: 


Be® + He* —> Cl? +7 (8.1) 


A source can be produced by mixing a natural a-emitter such as radium 

or polonium with beryllium powder, and in fact the Ra-a-Be source is 

widely used. The energy of the emitted neutron depends on (a) its 

direction with respect to the incoming particle, (6) the excitation energy 

of the residual carbon nucleus, and (c) the energy of the a-particle at 
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the instant of collision. As a result, the neutrons from such a source 
are distributed in energy between zero and about 13 Mev, with a broad 
maximum near 4 Mev (1). Ra-aBe sources yield about 10’ neutrons/sec 
per gram of radium. 

Neutrons can also be produced by the action of y-rays on deuterium 
or beryllium. Sources based on (vy, 7) reactions are called photoneutron 
sources. A Rary-Be source, consisting of a quantity of radium sur- 
rounded by beryllium metal has a neutron yield about one-tenth of 
that of a Ra-a-Be source and a spectrum extending up to about 1 Mev. 

Particle accelerators can be used to produce intense fluxes of neutrons. 
Beams of protons striking targets of tritium or lithium, and deuterons 
striking deuterium, tritium or beryllium all generate neutrons. Since 
all of these reactions can take place at less than 2 Mev (and some at 
much lower energies) Van de Graaff and Cockcroft-Walton accelerators 
can be used effectively as neutron generators. Both types of generator 
produce nearly mono-energetic beams of charged particles, and when 
these strike thin targets of any of the above-mentioned elements except 
beryllium, the neutrons emitted at any particular angle are also mono- 
energetic. (Beryllium is an exception because the residual nucleus can 
be left in an excited state, with a consequent reduction in the energy 
of the emitted neutron.) 

At much higher energies mono-energetic neutrons can be produced 
by the “stripping’’ process, which occurs when high-energy deuterons 
encounter matter. Some of the deuterons collide with nuclei of the 
target material in such a way that the proton is stopped while the neutron 
proceeds on its way, carrying about half the energy of the original 
deuteron. 

The most copious supply of neutrons is provided by the fission process 
as it occurs in a nuclear reactor. Fission neutrons are produced with a 
continuous distribution of energies with an average of about 2 Mev and 
an upper limit about 10 Mev. 

Flux measurements. Unlike charged particles, neutrons do not 
produce ionization in matter, and hence they cannot be detected directly 
in ionization chambers or counters. The measurement of neutron flux 
must therefore be based on interactions with nuclei which result in 
ionizing particles. The chief of these are nuclear reactions leading to 
the emission of a charged particle or to fission, capture leading to a 
radioactive nuclide which can then decay by the emission of ionizing 
radiation, and elastic scattering in which the recoil nucleus has enough 
energy to produce ionization. 

For any of the processes listed above the general principles of an 
absolute flux measurement are the same. Suppose a volume V of the 
material in which reactions are to be detected is placed in a flux distribu- 
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tion of neutrons ¢(v) (i.e., 6(v) dv is the flux of neutrons with velocities 
in dv at v). If the amount of material is small enough not to disturb the 
flux, the number of reactions of a particular kind produced per second 
will be 


R=NV f alew-ae (8.2) 


where N is the number of nuclei per unit volume and a,(v) is the cross 
section per nucleus at the velocity v for the reaction being observed. 
Since this general expression for the reaction rate involves an integral 
of the product of the flux and the reaction cross section, it does not 
give the total flux directly. However, in a number of special cases that 
occur frequently in practice the situation is simpler: 

(1) When the neutrons all have the same velocity, Eq. 8.2 reduces to 


R = NV¢o,(v) (8.3) 


If the cross section at the velocity v is known, a measurement of the 
rate at which the reactions occur yields the neutron flux. 

(2) At low energies the reaction cross section often varies inversely 
as the neutron velocity, i.e., 


Vo 
o,(v) = o,(v9) . (8.4) 
where oa,(vo) is the cross section at some particular velocity v9. Then 


R=NV f $(v)¢-(v9) = de 


and since ¢@ = nv, 


R = NVa,(v9)vo9 J n(v) du (8.5) 


where f n(v) dv is the total density of neutrons of all velocities. The 


measurement thus gives the total neutron density directly. If the 
spectral distribution is known—as it is, for example, for thermal neutrons 
—the total flux can be calculated. 

(3) In the case of a neutron beam striking an “opaque” detector, 
i.e., one in which all the neutrons are absorbed to produce the measured 
reaction, the situation is also simplified. If the area of the beam is A, 
then 


R=A f o(v) dv (8.6) 


and the total flux in the beam can be obtained without a knowledge of 
the cross section. 
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Thermal-neutron fluxes are more readily measured than those at 
higher energies because a number of elements have large reaction cross 
sections at thermal energies, and several of these cross sections are ac- 
curately know. Thus, for example, the reaction B® +n > Li? + a 
has a high 1/v cross section, and the reaction releases over 2 Mev in 
the kinetic energies of the a-particle and lithium nucleus. In a pro- 
portional counter filled with a known amount of BFs3 gas, each neutron 
reaction leads to a large pulse of ionization, and such a counter can be 
used to measure thermal-neutron fluxes. 

Thermal-neutron fluxes have also been measured with thin foils of 
indium or gold. These elements have large cross sections for the capture 
of a thermal neutron to form a radioactive isotope. The number of 
such radioactive nuclei produced in a given time can be determined by 
absolute counting techniques and the neutron flux deduced from the 
known capture cross section. 

Fast neutrons can be measured by first slowing them down to thermal 
energies in 8 moderator such as paraffin or water and then measuring 
the resulting thermal flux by one of the methods described above (2). 
Alternatively, the neutron-proton scattering process can be used. Here 
the number of scatterings with hydrogen nuclei in an appropriate 
material is measured by detecting the ionizing effects of the recoil 
protons, either in a counter (3, 4) or in a photographic emulsion. Not 
all the recoils can be detected, since their energies extend down to zero, 
but under controlled conditions the fraction missed can be estimated and 
taken into account. 

A more extensive discussion of neutron measurements can be found 
in @ review by Feld (5). 

Spectral measurements. Considerable work has been done on the 
measurement of slow-neutron spectra. Two common instruments for 
this purpose are the mechanical velocity selector and the crystal spec- 
trometer. In its simplest form the mechanical velocity selector or 
“chopper” consists of a pair of disks of neutron-absorbing material 
(usually cadmium) mounted at opposite ends of a rotating shaft. Each 
disk has a number of radial slits, those in the second disk being turned 
through a small angle with respect to those in the first. The neutron 
beam is directed parallel to the shaft at such a distance from the axis 
that it encounters the disks near their outer rims. A neutron passing 
through a slit in the first disk is stopped by the second disk unless its 
velocity is such that it arrives at the instant a slit is passing through the 
beam axis. The system thus transmits only neutrons in a narrow range 
of velocities determined by the speed of rotation of the disks and the 
angle between the sets of slits. By varying either of these parameters 
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and counting the transmitted neutrons, it is possible to measure the 
spectral distribution of a slow-neutron beam. A velocity selector using 
cadmium disks can only be used for energies up to the cadmium cut-off 
energy of about 0.4 ev. Other designs have extended this range beyond 
1000 ev. 

The crystal spectrometer makes use of the fact that a neutron of 
energy E has associated with it a wavelength \ = h/(2ME)”, where M 
is the neutron mass. In the case of slow neutrons this wavelength is of 
the order of the spacing between the planes of atoms in a crystal; hence 
slow neutrons undergo Bragg reflection in the same manner as x-rays, 
neutrons of only one energy being reflected for any particular angle be- 
tween the crystal plane and the beam direction. Varying the angle of 
the crystal makes possible an examination of the neutron flux as a 
function of energy. Crystal spectrometers are useful between about 
0.01 ev and 100 ev. 

The techniques outlined above cannot be applied to high-energy 
neutrons. Instead, use has been made of the recoil protons resulting 
from the scattering of the neutrons by the hydrogen in a photographic 
emulsion (6). If the track length of the recoil and the angle it makes 
with the direction of the neutron beam are measured, the energy of the 
incident neutron can be deduced. Measurements on a sufficient number 
of tracks—generally several hundred—then yield the distribution of 
neutron energies. 

Recoil-proton energies have also been measured in proportional- 
counter ‘“‘telescopes” in which the protons ejected in a particular direc- 
tion from a suitable radiator are detected by a pair of proportional 
counters operated in coincidence (7, 8). At energies below about 1 Mev 
the proton energy can be determined from the sum of the pulse sizes 
in the two proportional counters; at higher energies the pulse height 
produced in a scintillation counter has been used for the energy de- 
termination. 

Another instrument that shows promise for neutron spectrometry 
at high energies is the Li®I(Eu) scintillation counter (9, 10). Neutrons 
of all energies can initiate the reaction Li®(n, a)T, with an energy release 
in the scintillator of 4.78 Mev plus the neutron energy. Neutrons 
with energies above 1 Mev produce pulses that can readily be dis- 
tinguished from the slow-neutron pulses at 4.78 Mev, and since the 
reaction cross section is known as a function of energy, the pulse-height 
distribution from the counter can be used to deduce the fast-neutron 
spectrum. Although the resolution of this device is not great (about 
10% at 5 Mev), it has the the advantage of a relatively high detection 
efficiency that is independent of neutron direction. 
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Finally, the system of threshold detectors described on p. 104 can be 
used to obtain a rough estimate of the neutron spectrum. 

First-collision dose. Neutrons deposit energy in matter through 
the recoil nuclei from elastic scatterings, through charged particles 
emitted in nuclear reactions, and, to a degree depending on the dimen- 
sions of the medium, through y-rays following capture or inelastic 
scattering. Between thermal energies and about 10 Mev, by far the 
most important of these mechanisms is elastic scattering. At higher 
energies inelastic scattering and nuclear reactions become more im- 
portant, and at thermal energies the main contribution is from capture 
y-rays and, in some materials, the radiations from nuclei made radio- 
active by neutron capture. 

In general the relationship between absorbed dose and neutron flux 
is difficult to calculate, since y-rays produced at one point can contribute 
to the dose at other points where the neutron flux is quite different, and 
since scattering alters the energy spectrum and angular distribution of 
the neutrons. Some information can be gained from consideration of 
the limiting case where the neutrons interact with a layer of material 
so thin that more than one collision by a single neutron is unlikely. 
(The layer must, however, be thick compared to the ranges of the 
secondary charged particles produced by the neutrons, but this second 
condition can be readily met without violating the first.) The absorbed 
dose calculated for these conditions (i.e., the dose resulting from “first 
collisions”) is called the first-collision dose.* For neutrons of energy E 
it is given by 

D= 0.01¢tEp ZN orf: (8.7) 


where D is the absorbed dose in rads, ¢ is the neutron flux, ¢ is the time 
during which the material is exposed to this flux, p is the density, N; 
is the number of atoms of type k per cm?®, o; is the cross section for 
interaction with atoms of this type, and f; is the average fraction of the 
neutron energy transferred to the medium in such interactions. (In 
the case of a nuclear reaction in which the energy of the outgoing particle 
is greater than that of the incoming neutron, f; is greater than one.) 
The sum extends over all the types of atoms in the material. (If more 


* At the time of writing (1958) a somewhat different interpretation of first-collision 
dose appears to be replacing that given above. According to this interpretation, 
first-collision dose is defined simply as the kinetic energy of the ionizing particles set 
in motion per unit mass by the primary radiation. This definition is analogous to 
that of exposure dose, except that it refers to the energy of the secondary particles 
rather than their ionization. It would not alter Eq. 8.7 or Fig. 8.1. 
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than one kind of interaction can take place with the same type of atom, 
one may write 


Onfk = Djor sfr (8.8) 


where ;o; is the cross section for reactions of type j with atoms of type 
k, and ,f, is the average fraction of the neutron energy imparted to the 
medium in such a reaction.) 

The first-collision dose per neutron/cm? for tissue is shown as a func- 
tion of neutron energy in Fig. 8.1, which is based on the tissue composi- 
tion and values of o,f; given in reference 11. The figure also shows the 
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Fig. 8.1. First-collision curve for neutrons in tissue (solid curve) and total dose at 
surface of 30-cm slab (dashed curve). 


dose at the surface of a slab of tissue 30 cm thick, calculated to include 
the effects of multiple collisions in the medium (12). The same calcula- 
tions have been extended to various depths in the slab of tissue. 

Measurement of absorbed dose by cavity ionization. As in the 
case of x-rays, the ionization in a gas-filled cavity in a solid medium can 
be used to determine the absorbed dose due to neutrons interacting in 
the medium (13). When the cavity is small enough so that the secondary 
particles lose only small fractions of their energies in crossing it, the 
relationship between absorbed dose and ionization is given by the 
Bragg—Gray equation (Eq. 6.8). The stopping-power ratio appearing 
in the equation must be calculated for the actual spectrum of secondary 
particles crossing the cavity, mainly recoil protons in most cases of 
interest. Wg is the mean ionization potential for the particles in the 
gas, and may differ slightly from that for electrons (see Table 2.4). 
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The most common problem in neutron dosimetry is that of finding 
the absorbed dose at a point in tissue. If the wall material of the 
cavity differs in atomic composition from the surrounding tissue, the 
dose in the tissue is given by 


. (2N pone) 
(ZNicifi)u 


where Dy 1s the dose in the wall material, as given by Eq. 6.8. The 
bars over the quantities o,f; indicate averages over the neutron spectrum 
at the point in question. Unless the wall material is closely similar to 
tissue in composition—particularly in the amount of hydrogen it con- 
tains—the conversion factor differs considerably from unity. Since the 
cross sections and the neutron spectrum are often not known accurately, 
the value of the conversion factor tends to be uncertain when the wall 
differs greatly from tissue. Because of this difficulty, and because of 
the restriction imposed on the size of the cavity in order for the Bragg— 
Gray condition to be valid for walls and gas of different compositions, 
there is a considerable advantage in working with a tissue-equivalent 
gas enclosed in a tissue-equivalent wall. In this case the dose is given to 
a good degree of approximation by Eq. 6.9, and the gas volume may be 
as large as desired. The design and performance of tissue-equivalent 
chambers have been discussed in some detail by Rossi and Failla (14). 

One difficulty associated with the use of tissue-equivalent ionization 
chambers for neutron dosimetry is that such chambers respond to y- 
rays—in fact they measure the total absorbed dose due to both neutrons 
and y-rays. If the biological effect of one rad of y-rays were the same 
as that of one rad of neutron radiation, this might be a desirable property, 
but since the RBE of neutrons is several times greater than that of 
y-rays, the doses of the two types of radiation must be measured sepa- 
rately. Since neutrons are almost invariably accompanied by 7y-radia- 
tion, the problem is an integral part of neutron dosimetry. 

This difficulty can be at least partly overcome by the use of a second 
chamber of the same size and shape as the tissue-equivalent one but 
containing less hydrogen in the walls and gas. Such a chamber will 
have almost the same response to y-rays as the tissue-equivalent chamber 
but a reduced response to neutrons. For example, a graphite or teflon 
chamber filled with carbon dioxide has only about a tenth the response 
to fast neutrons of a tissue-equivalent chamber (15, 16), and from the 
readings of the two, the neutron and y-ray doses can be found separately. 

Cavity ionization is the basis for another unit (now outmoded) that 
has frequently been used in fast neutron dosimetry. This is the ‘n unit,” 
which is defined so that 1 n of neutron radiation produces the same 


T= (8.9) 
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ionization in a 100-r Victoreen condenser chamber as 1 r of x-rays. 
The Victoreen chamber is designed solely for use with x-rays, and two 
chambers that give the same response to a given x-ray intensity may 
give somewhat different responses to a given neutron flux. Further- 
more the absorbed dose in tissue corresponding to 1 n depends somewhat 
on the neutron energy, ranging from about 2 rads at high energies to 
something over 3 rads at energies below 1 Mev. Consequently when 
the absorbed dose in rads can be determined directly, it constitutes a 
more satisfactory measure of neutron radiation than the number of n 
units. 

Measurement of absorbed dose by proportional counters. The 
problem of neutron dosimetry in the presence of y-rays would be simpli- 
fied by a device that could be operated in such a way as to be insensitive 
to y-rays while still detecting neutrons. The proportional counter is 
such a device. Since the heavy charged particles set in motion by fast 
neutrons have much higher specific ionizations than the electrons set in 
motion by y-rays, they produce much larger pulses in traversing the 
gas of a proportional counter. If the bias level of the associated circuit 
is set so that only these larger pulses are accepted and counted, the 
instrument can be made to ignore y-rays. 

The problem with a proportional counter is to make the response 
proportional to the absorbed dose in tissue. One solution (17) has been 
simply to count the total number of pulses greater than the bias voltage, 
and to adjust the response to resemble the first-collision dose curve in 
tissue by means of various radiators within the counting tube. Calcula- 
tion shows that the recoil protons from a slab of paraffin provide a 
first approximation to this curve, but give insufficient counts at the 
low-energy and high-energy extremes. Interactions in methane gas in 
the counter can be used to enhance the response at low energies, and a 
second paraffin radiator covered with a thin aluminum absorber gives 
extra counts at high energies. The relative contributions from the dif- 
ferent radiators can be adjusted to make the response between 0.2 and 
10 Mev “‘tissue-equivalent” to within about 10%. A disadvantage of 
this system is that the sensitivity of the counter is directional, varying 
by as much as a factor of two depending on its orientation in a neutron 
beam. 

The other means of achieving tissue-equivalent response with a pro- 
portional counter has been to use a counter lined with polyethylene and 
filled with ethylene gas, and to integrate the pulse sizes so that, for 
example, a 2-volt pulse contributes twice as much to the reading of the 
instrument as a l-volt pulse (18, 19). The response of polyethylene to 
neutrons is similar to that of tissue, and the counter in effect makes use 
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of the Bragg-—Gray principle to give a reading proportional to the 
absorbed dose in polyethylene. Pulses below a certain size are again 
biased out to reduce the response to y-rays; this has the effect of dis- 
criminating also against low-energy neutrons, and the instrument is not 
reliable below about 0.2 Mev (16). 

It is also possible to operate a proportional counter in such a way that 
it detects y-rays and discriminates against neutrons (20). Such a 
counter is useful for measuring the y-radiation accompanying a neutron 
flux. 

Finally, the distribution of pulse heights from a spherical proportional 
counter operating at a very low pressure has been used to determine the 
distribution in specific ionization of the particles crossing the counter 
(21). This in turn provides a measure of absorbed dose as a function 
of LET, which is useful when the RBE depends strongly on LET. 

Absorbed dose from photographic emulsions. The absorbed 
dose due to fast neutrons can be determined by exposing a suitable 
photographic emulsion to the neutron flux and measuring the number 
and length of the proton recoil tracks produced. From a knowledge of 
the difference in composition between the emulsion and tissue, the 
absorbed dose in the latter can then be calculated. Such a procedure is 
obviously time-consuming. 

Threshold detectors. A series of threshold detectors can be used to 
determine the neutron flux in each of a series of energy intervals (22, 23), 
and from this information the first-collision dose can be calculated. 
The procedure is roughly as follows. First the total flux of all energies 
above about 1 kev is measured by observing the fission rate in a sample 
of Pu?*® shielded by B!°; Pu?*® shielded in this way has a nearly constant 
cross section down to 1 kev. Next the fission rate induced in Np?” is 
observed ; since this reaction hasa threshold at 0.7 Mev, it gives a measure 
of the flux above 0.7 Mev. The difference between this and the first 
measurement yields the flux between 1 kev and 0.7 Mev. Similarly 
measurements of the fission rate in U?*® (threshold 1.5 Mev) and the 
(n, p) reaction in S** (threshold 2.5 Mev) lead to the fluxes in the energy 
intervals 0.7 to 1.5 Mev, 1.5 to 2.5 Mev, and above 2.5 Mev. This 
method has the advantage that it is completely insensitive to y-rays. 


REFERENCES 


Teucher, Z. f. Phystk 126, 410 (1949). 
Hanson and J. L. McKibben, Phys. Rev. 72, 673 (1947). 


1. M. 
2. A. O. 
3. W. D. Allen and A. T. G. Ferguson, Proc. Phys. Soc. A70, 639 (1957). 


NEUTRON DOSIMETRY 105 


8. J. Bame, Jr., E. Haddad, J. E. Perry, Jr., and R. K. Smith, Rev. Sci. Instr. 28, 
997 (1957). 


. B. T. Feld in Experimental Nuclear Physics, edited by E. Segre, Wiley, New York, 


1953, Volume II, Part VII. 


. L. Rosen, Nucleonics 11, No. 7, 32 (1953); 11, No. 8, 38 (1953). 

. G. J. Perlow, Rev. Sci. Instr. 27, 460 (1956). 

. C. H. Johnson and C. C. Trail, ibid. 27, 468 (1956). 

. J. Schenck, Nature 171, 518 (1953). 

. R. B. Murray, Nuclear Instruments 2, 237 (1958). 

. Report of the International Commission on Radiological Units and Measure- 


ments (ICRU) 1956, Nat. Bur. Stand. Handbook 62 (1957). 


. W.S. Snyder and J. Neufeld, Brit. J. Radiol. 28, 342 (1955). 

. L. H. Gray, Proc. Camb. Phil. Soc. 40, 72 (1944). 

. H. H. Rossi and G. Failla, Nucleonics 14, No. 2, 32 (1956). 

. H. H. Rossi in Radiation Dosimetry, edited by Hine and Brownell, Academic 


Press, New York, 1956, Chapter 15. 


. H. H. Rossi, G. S. Hurst, W. A. Mills, and H. E. Hungerford, Nucleonics 13, 


No. 4, 46 (1955). 


. G.S. Hurst, R. H. Ritchie, and H. N. Wilson, Rev. Sct. Instr. 22, 981 (1951). 

. G. 8. Hurst, Brit. J. Radiol. 27, 353 (1954). 

. E. B. Wagner and G. 8S. Hurst, Rev. Sct. Instr. 29, 153 (1958). 

. J. A. Auxier, G. S. Hurst, and R. E. Zedler, Health Physics 1, 1 (1958). 

. H. H. Rossi and W. Rosenzweig, Radiology 64, 404 (1955). 

. G. 8. Hurst, J. A. Harter, P. N. Henley, W. A. Mills, M. Slater, and P. W. 


Reinhardt, Rev. Sct. Instr. 27, 153 (1956). 


. G. 8. Hurst, W. A. Mills, F. P. Conte, and A. C. Upton, Radiation Research 4, 


49 (1956). 


CHAPTER 


9 


Secondary 
instruments 


The term secondary instrument will be applied to any detector that 
does not yield an absolute measure of radiation intensity or dose. In 
principle any device that gives a response which varies with the amount 
of radiation falling on it can be used as a secondary instrument, but to 
be useful in practice the device should have a response proportional to 
the field strength for radiation of any particular energy. In this chapter 
some of the types of secondary instrument that have been found most 
useful in dosimetry will be discussed, and methods of calibration of 
these instruments will be outlined. 


CAVITY CHAMBERS 


Cavity ionization chambers provide one important means for the 
absolute measurement of exposure and absorbed dose, and the principal 
features governing the response of such chambers to radiations of various 
kinds are discussed in Chapters 5, 6, 7, and 8. However, when the 
design of the chamber, the materials of construction, and the associated 
circuitry are chosen with a view to ruggedness, portability, and ease of 
use, the resulting instrument is generally no longer suited for absolute 
measurements and must be considered as a secondary device. Such 
secondary chambers are widely used in x-ray dosimetry. 

In most secondary chambers intended for use with x-rays the walls 
are made roughly air-equivalent and the filling gas is air, so that the 
response is roughly proportional to the exposure dose. Perhaps the 
most commonly used detector of this type is the Victoreen condenser 
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r-meter, which consists of an air-filled nylon or bakelite chamber, coated 
on the inside with a conducting graphite mixture and fitted with an 
aluminum central electrode. The chamber proper is mounted on a 
cylindrical stem which acts as a capacitor. Ionization in the air in the 
chamber allows charge to leak off the capacitor, and the amount of this 
loss, as measured by the change in voltage across the capacitor, provides 
a measure of the radiation exposure. The wall materials of the chambers 
are not exactly air-equivalent, but the response at low energies is made 
more or less ‘“‘flat’’ (i.e., independent of energy) by adjusting the length 
of the aluminum collecting electrode and by the addition of small 
amounts of heavier elements to the conducting coating. (Heavier ele- 
ments contribute additional photoelectrons at low energies, which 
augment the ionization.) In this way chambers intended for use in 
the conventional x-ray region can be made energy-independent to within 
a few per cent between 80 and 250 kvp. The heavy elements tend to 
enhance the response to some extent over this whole energy range, and 
consequently at higher energies the response falls off somewhat before 
becoming energy-independent again in the straight Compton region. 
Thus a Victoreen chamber calibrated for conventional x-ray energies 
gives too low a reading on high-energy x-rays or on y-rays of greater 
than about 200 kev, and a separate calibration must be carried out for 
this energy region. 

Of the secondary instruments in common use for the measurement of 
x-rays, cavity chambers have the advantage that they can be made to 
give readings most nearly proportional to exposure or absorbed dose, 
and they are generally employed when accurate measurements of these 
quantities are required. When radiations of very low intensity must be 
measured, on the other hand, the more sensitive Geiger and scintillation 
counters have an advantage. 

For neutron measurements, cavity chambers have found less use as 
secondary detectors than proportional and scintillation counters. 


GEIGER COUNTERS 


In some respects the response of a Geiger counter to x-rays is similar 
to that of a cavity ionization chamber, in that it is due mainly to sec- 
ondary electrons generated in the walls and crossing the sensitive 
volume. But whereas in an ionization chamber each electron con- 
tributes in proportion to the energy it loses in crossing the cavity, in the 
Geiger counter the contributions of all electrons are equal, each electron 
giving rise to just one pulse. For this reason the energy dependence of 
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the response of a Geiger counter differs somewhat from that of a cavity 
chamber. 

The efficiencies of Geiger counters with various cathode materials 
as determined experimentally, are shown as a function of photon energy 
in Fig. 9.1 (1, 2,3). (The quantity plotted is the percentage of the 
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Fig. 9.1. Efficiency of a Geiger counter as a function of photon energy for various 
cathode materials. 


photons striking the counter that produce pulses.) For a copper cathode 
the efficiency between 0.3 and 2.5 Mev is very nearly proportional to 
the photon energy, and consequently a copper counter has a response 
proportional to the intensity of radiation in this energy region. Such a 
counter cannot therefore be expected to give readings proportional to 
the exposure dose at different radiation energies, and in fact the number 
of counts corresponding to an exposure of one roentgen increases grad- 
ually with increasing photon energy above about 0.5 Mev. However, 
as shown (4) in Fig. 9.2, the response between about 0.25 and 2.5 Mev 
lies within +15% of the mean value for this region, and for some 
purposes this is good enough. At lower energies the onset of the photo- 
electric effect causes the response to rise sharply, so that Geiger counters 
are of little use in the conventional x-ray region. Of course, when 
radiation of only a single energy (or a constant energy distribution) is 
involved, one can achieve accurate measurements by simply calibrating 
the counter specifically for this energy. 
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Geiger counters are not suitable for measuring absorbed dose in 
media such as tissue, partly because they tend to be too large, but mainly 
because they differ too much from tissue in energy response. 

The chief advantage of the Geiger counter over the cavity ionization 
chamber is the much greater sensitivity that can be achieved with 
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Fig. 9.2. Response of Geiger counters with various cathode materials to one roentgen 
of x-radiation. 


simple circuitry. This sensitivity is useful in weak fields, but it also 
means that the counter cannot be used at the higher intensities frequently 
encountered in radiation work. 


SCINTILLATION COUNTERS 


A scintillation counter detects radiation by converting some of the 
energy lost by charged particles in a suitable scintillating material to 
light. The total amount of light produced during a period of irradiation 
is referred to as the light output. In practice a fraction of this light is 
converted to an electric current at a light-sensitive surface and this 
current, amplified by a photomultiplier, is what is actually measured. 
The relationship between the measured current and the amount of 
energy deposited in the scintillator depends on so many factors that it 
is best determined empirically by exposing the scintillator to a known 
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radiation field and measuring the resulting current. That is, the system 
must be calibrated at at least one radiation energy, and hence can only 
be used as a secondary device. 

The scintillation counter can be made even more sensitive to x-rays 
than the Geiger counter. The energy response depends on the nature 
and size of the phosphor, but under the right circumstances it can be 
made proportional to either the exposure or the absorbed dose in tissue. 

Measurement of exposure dose. In inorganic crystals such as 
sodium iodide the light produced by a secondary electron is strictly 
proportional to the energy expended by it in the crystal (5). In 
organic scintillators such as anthracene the situation is less clear-cut, 
but there is evidence to show that when electrons are set in motion 
within the scintillator they produce amounts of light proportional to 
their energies (6, 7). Thus, to the extent that this proportionality holds, 
the light output of a scintillator is proportional to the energy locally 
absorbed from the radiation field. 

If radiation of plane intensity J, is incident normally on a crystal of 
area A and thickness d for a time ¢, the light produced at a depth z will 
be 

aL = kI,tAe ey, dz (9.1) 


where k is a constant of proportionality. The exponential factor takes 
into account the attenuation of the radiation in traversing the thickness 
x. The use of the energy absorption coefficient, uo, in the exponential is 
based on the assumption that radiation scattered in the crystal behaves 
much like unscattered radiation when it comes to further energy losses 
in the crystal; since the scattered radiation makes only a small contribu- 
tion to the total light output, this assumption is quite adequate for 
purposes of calculation. The total light output is obtained by inte- 
grating over 2Z, 

L = kI,tA(1 — ee?) (9.2) 


This expression takes a simpler form if the dimensions of the crystal 
are small. For if u.d<« 1, Eq. 9.2 reduces to 


L = kIptApgd = kIptV ue (9.3) 


where V is the volume of the crystal. Evidently in this case the response 
does not depend on the direction from which the radiation strikes the 
crystal, and when radiation comes from more than one direction the 
plane intensity in Eq. 9.3 is simply replaced by the spherical intensity, 
q,: 

L = kI,tVue (9.4) 


Ill 
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The exposure corresponding to a spherical intensity J, is simply 
proportional to the intensity and to the energy absorption coefficient 
for air (Eq. 5.2): 

R = K’It(ua)a (9.5) 


Thus the ratio of the light output of the crystal to the exposure is pro- 
portional to the ratio of the energy absorption coefficients for the radia- 
tion in the crystal and in air: 
L = (ua) crystal (9. 6) 
R (ua)A 

The closer the atomic number of the crystal material is to that of air, 
the less the ratio in Eq. 9.6 will depend on energy. Most of the in- 
organic scintillators contain elements with atomic numbers much greater 
than that of air, and consequently they have a response that varies 
greatly with energy. On the other hand, organic scintillators such as 
anthracene (which consists entirely of carbon and hydrogen) are “‘air- 
equivalent” over a considerable range of energies. These conclusions 
are illustrated in Fig. 9.3, which shows the relative light output per 
roentgen as a function of energy for sodium iodide and anthracene. (The 
curves have been normalized to unity at 1 Mev.) The curve for sodium 
iodide shows no flat portion; that for anthracene is flat between about 
0.2 Mev, where the photoelectric effect has become negligible, and 
about 5 Mev, where pair production becomes appreciable. The effect of 
finite crystal thickness is shown by the third curve, calculated from 
Eq. 9.2 for a thickness of 2cm. Except for energies below 30 kev, the 
effect of finite crystal size on the response is slight and actually tends 
to improve the energy independence. 

The fall-off in the response of anthracene at low energies can be 
largely overcome by adding materials of higher atomic number to organic 
scintillators (8, 9). 

The above analysis assumes that no appreciable fraction of the 
secondary-electron energy produced in the crystal escapes through the 
boundaries, since such losses are strongly energy-dependent. If the 
size of the crystal is comparable to (or less than) the ranges of the most 
energetic secondary electrons to be encountered, escape losses must be 
compensated by surrounding the crystal with an “equilibrium” layer 
of material of similar composition, so that the energy lost through 
electrons leaving the crystal is balanced by energy gained through 
electrons entering the crystal from outside. 

Measurement of absorbed dose. A small organic crystal can also 
be used to measure relative values of absorbed dose. An analysis 
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similar to that leading to Eq. 9.6 gives for the ratio between the light 
output and the absorbed dose in a medium, M, 


L ee (Ha/ p) crystal 
D (Ha/ p)M 


The flattest energy response is obtained with a crystal approximating 
the medium of interest in composition. Thus unmodified anthracene 
gives a light output closely proportional to the absorbed dose in tissue 
between 0.2 and 3 Mev. 

Similarly an organic crystal can be used to measure the absorbed dose 
due to a beam of high-energy electrons. In this case the crystal must 
be small compared to the electron range. The ratio between the light 
output of the crystal and the absorbed dose in tissue is proportional to 
the ratio of the stopping powers; since this ratio changes very little 
with energy, the light output is closely proportional to the absorbed 
dose in tissue over a wide range of energies. 

Measurement of intensity. If, instead of being small, the crystal 
is so large that it absorbs practically all the energy falling on it, the 
light output is proportional to the intensity. For this purpose some of 
the inorganic crystals are superior to organic crystals, since they have 
higher densities and contain elements of higher atomic number. For 
example, a 5-cm cube of sodium iodide exposed to 200-kev photons will 
absorb all but a few per cent of the energy falling on it, and larger 
crystals are effective up to much higher energies. 

The intensity in a charged-particle beam can also be measured in this 
way. 

Pulse counting. The measurement of the integrated photomultiplier 
current due to the light being produced in the crystal is only one way 
of using scintillation detectors. Another is to count the individual 
pulses of current arising from the interactions of individual photons in 
the crystal. The number of such pulses is proportional to the number of 
photons interacting in the crystal, rather than to the energy absorbed. 

The total count, C, due to radiation of intensity J, and photon energy 
hy falling on a small crystal of volume V for a time ¢ is 

ItV 


C = k’ —uzu (9.8) 
hy 


(9.7) 


where pz is the total absorption coefficient at the energy hy. The factor 
hy in the denominator rules out the possibility of any proportionality 
between count and exposure or absorbed dose, even over a restricted 
range of energies. 
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If the crystal is large enough so that practically every photon incident 
upon it interacts to give a pulse, the observed counting rate is propor- 
tional to the total number of photons arriving per second, i.e., to the 
photon flux. If the pulse-height distribution is measured as well as 
the total count, the spectral distribution of the radiation can be de- 
termined (see Chapter 3), and this makes it possible to calculate the 
intensity from the photon flux (10). 

Charged-particle flux can also be determined from the counting rate in 
a crystal. Since a charged particle loses some energy even in a thin 
layer of matter, the crystal can be small and still give a pulse of light 
for each particle striking it. 

Neutron dosimetry. In principle an organic crystal should provide 
an excellent means of measuring fast-neutron dose in tissue, since the 
light output from the crystal and the absorbed dose in tissue are both 
due mainly to recoil protons. In practice, however, the flux of y-rays 
that generally accompanies a neutron flux greatly reduces the usefulness 
of this method, since y-rays make a much greater contribution to the 
light output in relation to their biological effect than do neutrons. A 
number of ingenious methods of discriminating against y-rays while 
still retaining at least a rough tissue equivalence for neutrons have been 
devised (11). For example, zinc sulphide powder, which is a very ef- 
ficient scintillator, can be mixed with a substance with a high hydrogen 
content, such as Lucite or paraffin (12). Some of the recoil protons 
from the hydrogenous material pass through grains of zinc sulphide and 
produce light pulses, and because the specific energy loss of a recoil 
proton is much greater than that of a recoil electron from a y-ray inter- 
action, these pulses are much larger than those due to y-rays. By 
counting only pulses greater than a certain size, it is thus possible to 
discriminate against y-rays and detect only neutrons. Unfortunately 
the response of such a system is not tissue-equivalent, and the detector 
must be calibrated for the particular neutron spectrum on which it is 
to be used. 


PHOTOGRAPHIC FILMS 


A charged particle passing through a grain of silver bromide in a 
photographic emulsion will generally cause changes which will result 
in the conversion of the grain to atomic silver when the film is developed. 
Such developed grains cause a blackening of the emulsion which can 
be measured as an increase in photographic density. In this way it is 
possible to detect electrons and heavy particles and, through the effects 
of their secondary electrons, x-rays. 
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The extent of the blackening for a given exposure depends on, among 
other things, the type and thickness of the emulsion, the nature, energy 
and angular distribution of the radiation, and the conditions of develop- 
ment. Although no precise relationship between film density and condi- 
tions of irradiation has been worked out, a few regularities may be 
mentioned. For small exposures the density in many types of emulsion 
is proportional to the energy absorbed from any particular type and 
energy of radiation; a different type or energy of radiation may have a 
different constant of proportionality. However, since the maximum 
density obtainable is limited by the amount of silver bromide in the 
emulsion, this proportionality necessarily breaks down as the dose in- 
creases. For high-energy §-rays the response does not depend greatly 
on energy; however, it falls off badly at low energies. The response per 
roentgen of x-rays is virtually independent of photon energy above 0.3 
Mev, but below this energy it rises several fold as photoelectric inter- 
actions in the silver come into play. A typical response curve for x-rays 
is shown in Fig. 9.4 (13, 14), where the ordinate is proportional to the 
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Fig. 9.4. Typical curve of film response to x-rays. The ordinate is proportional to 
the reciprocal of the number of roentgens required to produce a standard density. 


reciprocal of the number of roentgens required to produce a particular 
density. The response at low energies can be flattened considerably by 
the addition of filters which selectively absorb the low-energy com- 
ponents (13), although this procedure tends to make the response more 
directional. 

Exposures of 0.1 to 1000 r can be measured with suitable films— 
roughly the same range as covered by condenser ionization chambers. 
Films are used to advantage for dosimetry where a very thin detector 
is desirable in order not to disturb the experimental conditions, or where 
a compact and rugged detector is required. 
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CHEMICAL DOSIMETERS 


At the opposite end of the sensitivity scale to Geiger and scintillation 
counters are the chemical dosimeters. Many substances undergo 
chemical changes under the action of ionizing radiations; this is par- 
ticularly true of chemicals in aqueous solution, where the decomposi- 
tion of the water itself contributes to the reaction. If the extent of 
the chemical change can be conveniently measured, the reaction can be 
used for dosimetry. 

A widely used reaction will serve as an example of the properties of 
chemical dosimeters: the oxidation of ferrous sulphate to ferric sulphate 
(15). The concentration of ferric ions produced by irradiation can be 
determined by measuring the absorption of light of a particular wave- 
length in a standard thickness of solution with the aid of a spectro- 
photometer. The yield of ferric ions is proportional to the absorbed 
dose for doses between about 10? and 10° rads. The yield is substantially 
independent of radiation energy for x-rays above about 200 kev., but 
tends to fall off at higher LET values. Thus for measurements with 
x-rays and y-rays the ferrous-ferric system need not be calibrated for 
each individual application; instead, the experimentally determined 
yield of 15.6 + 0.2 ions per 100 ev expended in the solution (16) may 
be used over a wide range of energies to deduce the absorbed dose in the 
solution. Since the composition of an aqueous solution is similar to 
that of tissue, a system of this kind can be used to give a direct measure 
of absorbed dose in tissue over a wide range of radiation energies. 

The low sensitivity of the ferrous-ferric system restricts its use to the 
measurement of large doses. However, a number of other systems have 
been investigated, some of which can measure doses down to a few 
rads. The principal features of a number of chemical systems are 
described in a review by Taplin (16). 


METHODS OF CALIBRATION 


Secondary instruments must be calibrated against a primary measur- 
ing system, either by direct comparison or through the medium of a 
secondary standard. The method of calibration depends, of course, on 
the type and energy of radiation to be measured. 

X-rays. When the greatest possible accuracy is required, instruments 
to be used for the measurement of x-ray exposures at voltages up to 
300 kvp are generally calibrated directly against a free-air chamber at 
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one of the standardizing laboratories. At these energies few secondary 
instruments can be relied upon to be energy-independent, and calibra- 
tions must be made at a number of tube voltages and filtrations covering 
the range of qualities over which the instrument is to be used. An 
instrument so calibrated can then be used in turn to calibrate other 
instruments to somewhat less accuracy. An instrument chosen for the 
constancy of its response * over extended periods of time and used in 
this way to calibrate other instruments is sometimes referred to as a 
secondary standard. 

As pointed out in Chapter 5, existing free-air chambers are not 
suitable for calibrations of x-rays much above 300 kvp. At higher 
energies such as those encountered in y-rays, comparison must be made 
with a cavity chamber which has been designed for absolute measure- 
ments. If the secondary instrument is itself a form of cavity chamber 
and is to be used for measurements on 7-rays, a single calibration at a 
single convenient energy will usually suffice, since any chamber con- 
structed of low-Z materials is practically air-equivalent over the normal 
range of y-ray energies. In the case of other devices such as Geiger 
and scintillation counters the desirability of calibrating at more than 
one energy will depend on the accuracy required. 

Calibrations at y-ray energies are often made with the aid of a radio- 
active source with a known y-ray output, the output being based in turn 
on an absolute cavity-chamber measurement. One of the sources most 
commonly used for this purpose is radium. The output of a milligram 
of radium in equilibrium with its decay products and encapsulated in 
0.5 mm of platinum, has been determined in a number of careful experi- 
ments (18): the value adopted by the national laboratories in North 
America is 8.25 r/hour at 1 cm (19). A radium source encapsulated in 
this way and containing an accurately known amount of radium thus 
provides a known 7-ray field in which secondary instruments can be 
calibrated. If the wall thickness of the capsule differs from the standard 
0.5 mm of platinum, the output can still be calculated by taking into 
account the difference in absorption between the actual capsule and the 
standard one (20). 

One advantage of radium as a secondary standard source is that the 
half-life is so long (1600 years) that corrections for decay need be made 
only at intervals of several years. By contrast, cobalt-60, which is also 
sometimes used as a secondary standard, has a half-life of 5.2 years, 
and hence decays at the rate of about 1% per month. 

* The response of a detector can be checked for constancy by exposing it periodi- 


cally, under fixed geometrical conditions, to a source of constant radiation output, 
such as radium. 
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Many artificially produced radionuclides have simple decay schemes 
from which the y-ray output per millicurie can be calculated. However, 
such figures are of limited value for sources of the size used in roentgen 
calibrations, since it is generally more difficult to determine the activity 
in millicuries than to measure the output directly. 

Beta-rays. Because §-rays are so strongly absorbed in matter, most 
secondary detectors can give little more than an order of magnitude for 
the B-ray dose. Detectors such as thin-walled ionization chambers or 
thinly wrapped photographic films can be calibrated for f-rays of a 
particular energy by placing them at a fixed point in relation to a f-ray 
source that has previously been calibrated with an absolute device such 
as an extrapolation chamber. 

A rough calibration of a thin-walled ionization chamber for high- 
energy §-rays can be achieved by adding an equilibrium wall thickness 
and exposing the chamber in a known 7-ray field. The number of rads 
in air corresponding to the observed ionization can be calculated from 
the parameters of the y-ray field, and the resulting conversion factor can 
then be applied to §-ray measurements made with the thin wall. If 
a significant fraction of the §-ray spectrum to be measured consists of 
electrons too slow to penetrate the chamber wall, the reading obtained 
will be too low. 

Neutrons. Secondary neutron detectors are often calibrated in a 
standard neutron flux. Standard thermal-neutron fluxes have been 
established by placing one or more neutron sources in a moderating 
medium and calibrating the resulting thermal flux at a standard position; 
the standardization is carried out by observing the rate at which nuclear 
reactions occur in an element for which the thermal-neutron reaction 
cross section is accurately known (21). A known flux of fast neutrons 
can be obtained from a neutron source for which the total neutron 
emission has been determined by one of the standard methods (22, 23). 
Such sources have complex energy spectra, and the absorbed dose in 
tissue must be calculated on the basis of this spectrum before the flux 
can be used for the calibration of fast-neutron dosimeters. 
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Cavity ionization chamber, comparison 
against calorimeter, 79 
equilibrium wall thickness of, 63, 64 
exposure dose by, 61-70 
for neutrons, 101-103 
stopping-power correction for, 65-69 
Cerenkov radiation, 9 
Charged particles, dosimetry of, 86-93 
interactions of, 8-26 
Chemical dosimeters, 116 
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Elastic scattering, of electrons, 9, 18-20 
of neutrons, 40, 41 
Electrons, absorbed dose from, 88-90 
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energy loss by radiation, 16-18 
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Energy loss by radiation 
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ments, 118 
for electron measurements, 87, 88 
for heavy-particle measurements, 93 
tissue-equivalent, 77, 102 
(See also Cavity ionization chamber, 
Free-air chamber) 


Klein-Nishina formula, 29, 68 


LET, 14, 92 
Linear accelerator, 86, 92 
Linear energy transfer, 14, 92 


Mean excitation potential, 10, 11, 12 
for air, 11 
for solids, 11 
table of values, 12 
Multiple scattering, 19 
values for electrons, 19 


n unit, 102, 103 
Neutron density, 6 
Neutrons, capture, 42 
dosimetry, 95-104 
elastic scattering, 40, 41 
flux measurements, 96-98 
inelastic scattering, 41 
interactions with matter, 40-42 
nuclear reactions, 41 
secondary detectors for, 118 
slowing down, 42 
sources, 95, 96 
standard flux, 118 
Nonparticipating electrons, 12 
Number-distance curve, 20, 22 


Pair production, 26, 33-34 
Penetration of charged particles, 20-24 
Photoelectric effect, 26-28 
Photographic emulsions, as secondary 
detectors, 114, 115 

for measuring neutron dose, 104 

for measuring neutron flux, 98 

for measuring neutron spectra, 99 

response per roentgen, 115 
Photoneutron sources, 96 
Photonuclear reactions, 34 
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Plane intensity, 3, 4 
Polarization effect (see Density effect) 
Positrons, annihilation of, 33 
production by x-rays, 33 
Proportional counters, for measuring 
neutron dose, 103, 104 
for measuring neutron flux, 98 
for measuring neutron spectra, 99 


Quantity of radiation, definition of, 4 
relationship to exposure dose, 56, 57 


Rad, 6, 56 
definition of, 6 
(See also Absorbed dose) 
Radiation field, 1-6 
Radiation loss, 16-18 
Radium, as standard y-ray source, 117 
y-ray output of, 117 
Range, extrapolated, 21 
maximum, 23 
mean, 21 
of charged particles, 20-24 
RBE, definition of, 92 
of y-rays, 92 
of heavy charged particles, 92 
of neutrons, 92, 95 
of x-rays, 92 
RBE dose, 92 
Relative biological effectiveness (see 
RBE) 
Relativistic increase in stopping power, 
9, 14 
Rem, definition of, 92 
Rep, 56 
Roentgen, definition of, 55 
(See also Exposure dose) 


Scattering, of electrons, 18-20, 22 
of neutrons, 40, 41 
of x-rays (see Compton effect) 
Scintillation counters, 2, 109-114 
measurement of absorbed dose by, 
112, 113 
measurement of exposure dose by, 
110-112 
measurement of intensity by, 113 
measurement of photon flux by, 114 
neutron dosimetry by, 114 
neutron spectrometry by, 99 
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Scintillation counters, response per roent- 
gen, 111, 112 
Scintillation spectrometry, of neutrons, 
99 
of x-rays, 47 
Secondary instruments, 106-118 
Sodium iodide, as scintillation detector, 
111, 112 
Specific energy loss, 9 
Spectral distribution, definition of, 1, 2 
measurement for neutrons, 98, 99 
measurement for x-rays, 44-48 
Spherical intensity, 2, 3, 4 
Stopping power, 9, 10, 11, 14 
for electrons, 10 
for heavy particles, 10 
mass, 11 
per electron, 11 
variation of with energy, 14 
Stopping-power ratio, 65-70, 75, 89, 
101 
Stripping reaction, 96 
Synchrotron, 86 


Thermistor, 50, 73 
Threshold detectors, 100, 104 
Tissue-equivalent ionization chamber, for 
neutron measurements, 102 
for x-ray measurements, 77 


Van de Graaff generator, 86, 92, 96 

Velocity selector (for neutrons), 98 

Victoreen condenser chamber, 106, 107 
neutron measurements with, 102, 103 
response to x-rays and 7-rays, 107 


W, 26, 76 
definition of, 26 
table of values, 26 


X-rays, characteristic, 45 
continuous spectrum, 16, 44 
determination of spectrum, 44—48 
measurement of absorbed dose, 72-84 
measurement of exposure dose, 55-70 
measurement of intensity, 50-54 
secondary detectors for, 116, 117 
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